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HECTIC”"  A  LOR  AN  Grid  Computation 
Program  for  Irregular  Terrain 


1.  INTRODUCTION 

Program  HECTIC  is  a  software  system  for  computing  time-difference  coor¬ 
dinates  (TDl  and  TD2)  produced  at  specified  geographic  coordinates  by  a  low  fre¬ 
quency  LORAN  triad.  The  program,  an  updated  version  of  HUFLOC,  ^  utilizes 

2 

algorithms  originally  developed  by  Johler  and  Berry  and  later  implemented  by 

3 

Horowitz  in  LORANCO  to  compute  secondary  phase  correction  factors  via  Hufford's 

4 

integral  equation.  A  data  base  is  required  which  specifies  the  earth  impedance 
and  terrain  elevation  at  points  within  the  geographical  region  of  interest.  As 
implemented  on  a  CDC  6600  computer,  HECTIC  incorporates  the  following  modifica¬ 
tions  to  the  LORANCO  system  configuration: 


(Received  for  publication  29  January  1980) 

1.  HUFLOC  (1977)  Computer  Program  from  USAF,  Tactical  LORAN  SPO,  \ 

Hanscom  AFB,  MA  01731.  | 

2.  Johler,  J.  R.  ,  and  Berry,  L.  A.  (1967)  Loran-D  Phase  Corrections  Over  I 

Inhomogeneous,  Irregular  Terrain,  ESSA  Technical  Report  lEK  59-ITSA  56,  | 

(Suptd.  of  Documents,  U.S.  Govt  Printing  Office,  Washington,  DC  20402).  | 

3.  Horowitz,  S.  (1977)  User's  Guide  for  ESD  LORAN  Grid  Prediction  Program,  | 

HADC-TR-77-407r"  | 

4.  Hufford,  G.  A.  (1952)  An  integral  equation  approach  to  the  problem  of  wave  | 

propagation  over  an  irregular  surface.  Quart.  Appl.  Math.  0(No.  4):391.  I 
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(1)  Output  has  beon  limitPd  to  data  npcessary  Cor  comparing  tiniP-difOrpncp 
predictions  with  nipasurpiiients, 

(2)  New  error  messages  and  error  recoveries  have  been  added, 

(3)  Database  handler  code  no  longer  is  hard-wired  Cor  a  specific  geographic 
area, 

(4)  'i'argpts  below  the  earth's  surface  are  rejected, 

(5)  The  input  card  stream  was  revised,  and 

(6)  The  use  of  blank  and  labelled  COMMON  was  made  more  economical. 


2.  COMPUTER  CODE 

2.1  Tiine-of-arrival  Function  TOA 

LORAN  time-difference  coordinates  TDl  and  TD2  for  a  prediction  grid  serviced 
by  a  3 -station  chain  are  given  by 

TDl  =  TOAl  -  TOAM 
TD2  =  TOA2  -  TOAM 

where  TOAl  and  TOA2  and  TOAM  are  the  signal  time-of-arrival  for  Slave  1, 

Slave  2,  and  Master  stations,  respectively.  Time-of-arrival  is  determined  in  the 
main  program  HECTIC  via 

10®  0 

TOA  =  5Ji  + - £.  +  eD  (ii  sec) 

c  w 

whe  re 

T]  =  atmospheric  index  of  refraction  at  the  surface  (=  1.  000338), 

X  =  geodesic  distance  from  transmitter  to  target's  projection  onto 

a  spheroidal  earth  of  equatorial  radius and  polar  radius  b^  (m), 
c  =  vacuum  velocity  of  light  (=  2.  997925  X  10®  m/  fxsec), 
u)  =  radian  frequency  (=  2nf  =  27r  X  10®  Hz), 

ED  =  emission  delay  ( 4  sec), 

c-nd 

<3^  =  secondary  phase  correction  factor, 

=  -  argd-'W)  +  k(r^  -  x). 


C^. 
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where 


!•'  =  induction  field  function  KZ  or  FH, 

VV  =  secondary  wave  attenuation  factor, 
k  =  wave  number  in  air  (ourj/c  m  ), 

and 

r  =  chordal  distance  (m)  from  transmitter  to  target's  projection  onto  a 

o  6 

spherical  earth  of  effective  radius  a/cv  where  a  =  6.36739  X  10  (m) 
and  (V  is  the  atmospheric  vertical  lapse  factor. 


2.2  Induction  Field  Function  F 

The  induction  field  is  formulated  in  subroutine  INDF  as  either 


a  a  sin  6 
FZ  ^ ^ -  + 


r  1 

■  3a  a  sin^  6 

ikD  ^  ( ikD  1 

FH  = 


a  sin  6 
r 

~2r-D~- 

o 


1  +  T 


ikD 


where 


a  +  h  , 
r 


h  =  target  altitude  (m). 


2  2 

=  a  +  a  -  2a  a  cos  6 , 
r  r 

=  x/a. 

=  ^  -  1,  and 

=  J  M  U  -  377  . 

>  o  o 


FZ  is  used  when  the  phase  of  the  vertical  electric  field  is  to  be  found,  and  FH  is 
used  when  the  phase  of  the  azimuthal  magnetic  field  is  required. 


Although  X  is  the  geodesic  distance  on  a  spheroidal  earth,  is  calculated  as  if 
X  we-re  the  distance  on  a  spherical  earth. 
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2.3  Seft>ndary  Wave  Atteimatic)ii  Kuiu-lioii  W’,  Surfaei’  <  ji«e 

\\  hen  the  target  is  on  the  eui-th’s  surface,  the  quantity  U  is  calculated  in  sub¬ 
routine  INFQ2F  from  the  following  integral  equation; 


\V(x)  =  1  -  ^  J  ds  •  VV(s)  ^ 


X 

s  (x-s) 


1  + 


3h 

s 


■1/2 


1  +  T 


ik  r„ 


‘2 

BIT 


s(x  -  s) 
'^l  ^2 


o 

X 


1  + 


9  h 

s 

FT-, 


1/2 


,-i k  r 


where 


s 


5 


s 


k/k2 
^  2 


9h 


s 

FT“ 


r 


r 


2 

1 


2 

r 

o 


=  distance  along  the  transmitter-to-target  geodesic  {0  -  s  -  x). 
=  normalized  impedance  at  point  s, 

=  (k/k2)<l  -  k^/k2^)^^^. 

=  relative  dielectric  constant, 

=  8.  854  X  10"^^  F/m, 

=  conductivity  (mho/m). 


=  terrain  slope  at  s, 

=  r  +  r„  -  r  , 

1  2  O 

=  a  +  a  -  2a  a  cos  (s/a), 
s  s 

=  a  -f  h  , 
s 

=  terrain  elevation  at  point  s  relative  to  terrain  elevation 

at  the  transmitter, 

2  2  .X  -  s  V 

=  a  +  a  -  2a  a  cos  ( - ), 

X  s  X  s  a 

=  a  +  h  , 


X 

^  terrain  elevation  at  point  x  relative  to  terrain  elevation 

at  the  transmitter, 

2  2  „  ,x. 

"  a  f  a  -  2a  a  cos  (— ), 

X  X  a 
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normal  derivative  of  r^  with  respect  to  the  earth's  surface. 


2.4  >e(oii<lary  Wave  .Vttrniialioii  hiiiitlioii  Klrvalnl  Caisr 

If  the  target  is  above  the  earth's  surface,  then  subroutine  also  cal¬ 

culates 


whore  W  (s)  is  the  previously  obtained  surface  function  W(x),  and  the  quantity  a 
which  enters  into  calculating  r  .  r  .  r  ,  r  and  3  r„/a  is  replaced  bv  a  =  a  +  h 

1^0  -ifl  "'J*  p* 

where  is  the  tai’get  altitude  relative  to  terrain  elevation  at  the  transmitter. 

:i.  SY.STKM  DKSKliN 

Program  HECTIC:  is  divided  into  3  modules;  Driver,  Database,  Handler  and 
■Secondary  Phase  Calculator. 

3.1  Driver  Module 

The  Driver  routine  Hl'.C  IIC  reads  data  from  input  cards,  checks  their  validity 
and  converts  them  from  user  to  internal  format  via  routines  C70KD.MS  and  COHKAD. 
It  .sets  up  the  three  geodesic  propagation  paths  from  the  stations  to  each  target  via 
routines  CIEDDI  and  CEDPT.S.  It  retrieves  impedance  and  terrain  data  from  the 
database  via  the  Database  Handler  routine  CiEOUEX.  It  also  calls  the  Secondary 
Phase  Calculator  routine  liMX^2E  and  prints  the  final  results  for  TDl  and  TD2. 


1  1 


is. 2  Dululiase  Maiidlfr  Mocl(ili> 

The  Database  Handler  routine  CiKOHKT  reads  the  database  tape  and  transfers 
the  impedanee  and  terrain  data  to  mass  storage  via  routine  C;1CTK1.\'.  It  also 
retrieves  these  data  fiom  mass  storage  via  loutines  CIDTDIA  ,  INDKX,  I\T,  and 
L  NPAtiK. 

iS.il  .''(■coiular)  lliuse  (Uileulaloi'  Module 

The  secondary  Phase  Calculator  routine  1NEQ2E  performs  initial  flat-earth 
calculations  via  routines  Fl.EAP  and  WEItP.  It  interpolates  the  impedance  and 
terrain  data  via  routines  C'NEVKKN,  GROUND,  INT,  and  OMCC^S.  Then  it  evaluates 
Ilufford's  integral  equation  and  computes  the  secondary  phase  via  routine  GANG 
and  INDP.  These  operations  are  diagrammed  in  the  System  Flow  Chart  inF'igure  1. 


4.  USKR’SGUIDK 


Program  flKCTlC  reads  input  cards  and  a  database  tape.  It  produces  printer 
output.  These  data  sets  are  described  below. 

4.1  Input  Curds 

The  input  file  consists  of  Control,  Map,  Station,  and  Target  cards.  Figure  2 
shows  a  sample  input  card  sequence. 

4.1.1  CONTROL  CARDS 

Cards  numbered  1  through  11  in  Table  1  are  program  control  cards.  These 
cards  specify  f,  tj  ,  a,  a  ,  b^,  c  and  a  which  were  defined  in  Section  2.  In  addi¬ 
tion,  the  flat-earth  distance  F^  1000  m)  specifies  how  far  along  the  path,  flat 
earth  theory  is  to  be  used  to  initialize  W(s)  before  round  earth  code  begins.  There 
are  also  two  control  switches:  one  selects  either  the  electric  field  or  the  magnetic 
field;  the  other  allows  the  selection  of  either  a  rough  or  a  smooth,  inhomogeneous 
earth  model.  (In  the  smooth  model  terrain  variations  are  ignored.  ) 

4.1.2  MAP  CARDS 

The  database  tape  format  does  not  include  geographic  coordinates  because 
HECTIC  automatically  keys  location  to  data  via  the  database  sequence  structure  (to 
be  described  in  Section  4.  2).  In  order  to  implement  this  keying  operation,  the 
user  supplies  the  location  of  the  first  data  point  via  Card  12,  and  the  location 
30”  east  and  30”  north  of  the  last  data  point  via  Card  13.  These  locations  are 
referred  to  as  the  southwest  and  northeast  map  corners,  respectively.  The  input 
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card  t'orniat  for  all  geographic  positions  is  ddd-nim -ss.  sss  A,  where  ddd  =  0  to  90° 
latitude,  or  0  to  180°  longitude,  mm  =  0  to  59  min,  ss.  sss  =  00.000  to  59.999  sec, 
and  A  =  .\,S,F  or  W.  The  map  corners  must  lie  on  integral  latitude  and  longitude 
lines: 

dd-00-00.  000  .V  or  S 
ddd-00-00. 000  F  or  W 

with  the  following  restrictions:  (1)  the  northeast  corner  cannot  lie  north  of  H9’N, 
(2)  the  southwest  corner  cannot  lie  south  of  89°S,  (3)  the  niap  cannot  straiidle  the 
180°  meridian,  and  (4)  llie  minimum  map  region  is  a  zone  1°  latitude  by  2’  longi¬ 
tude,  with  larger  n  ap  regions  being  multiples  of  this  basic  zone. 


CALCULATOR 


I'igure  1. 


•Svsten.  I'low  t'iiar' 
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Figure  2.  Sample  Input  Card  Images 


Table  1.  Input  Card  Formats 


Card 

Card 

Read 

Number 

Column(s) 

Format 

Data  Entry 

1 

5 

15 

1 

6-16 

ElO.  0 

Loran  transmission  frequency  f  (kHz) 

2 

5 

15 

2 

6-16 

ElO.  0 

Air  index  of  refraction  at  ground  t] 

3 

5 

15 

3 

6-16 

ElO.  0 

Spherical  earth  radius  a  (m) 

4 

5 

15 

4 

6-16 

ElO.  0 

Spheroidal  earth  equatorial  radius  (m) 

5 

5 

15 

5 

6-16 

ElO.  0 

Spheroidal  earth  polar  radius  (m) 

6 

5 

15 

9 

6-16 

ElO.  0 

Vacuum  speed  of  light  c  (km//zsec)  j 

7 

4-5 

15 

11  j 

6-16 

ElO.  0 

Flat-earth  distance  F  (m) 

8 

4-5 

15 

\2 

6-16 

ElO.  0 

Vertical  lapse  factor  a 

9 

4-5 

15 

14 

6-  16 

ElO.  0 

1.  0  for  electric  field  prediction; 

0.  0  for  magnetic  field  prediction. 
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I'ublp  1.  Input  Oai  ii  l''ornial.s  (t.'onl) 


Card 

t.'ard 

Head 

.Number 

( 'olumrds ) 

Format 

Data  Fntrv 

10 

4-5 

15 

15 

6-16 

FIO.  0 

0.0  for  data  base  rough  earth  terrain; 

1.0  for  smooth  earth. 

1 1 

1-5 

lllank.  Terminates  t'ontrol  Card 
stream. 

12 

2-12 

AlO 

MAP  SW. 

13-17 

15 

Map  .SW  corner  latitude  degrees. 

18-20 

13 

Map  .SW  corner  latitude  minutes. 

21-27 

K7.  3 

MAP  SW  corner  latitude  seconds. 

28 

A1 

Map  SW  corner  hemisphere,  N  or  S. 

29-33 

15 

Map  SW  corner  longitude  degrees. 

34-36 

13 

Map  SW  corner  longitude  minutes. 

37  -43 

F7.  3 

Map  SW  corner  longitude  seconds. 

44 

A1 

Map  SW  corner  hemisphere,  F  or  W. 

13 

2-12 

AlO 

MAP  NE. 

13-44 

Map  NE  corner.  Use  the  formats  of 

Card  12. 

14 

2-12 

AlO 

MASTER. 

13-44 

Master  Station  position.  Use  the 
formats  of  Card  12. 

15 

2-12 

AlO 

SLAVE  1. 

13-44 

Slave  1  position.  Use  the  formats 
of  Card  12. 

45  -53 

F9.  2 

Slave  1  emission  delav  (/jsec). 

16 

2-12 

AlO 

SLAVE  2. 

13-44 

Slave  2  position.  Use  the  formats  of 

Card  12. 

45-53 

F9.  2 

Slave  2  emission  delav  (/isec) 

17 

2-12 

AlO 

Target  identifier. 

13-44 

Target  position.  Use  the  formats 
of  Card  12. 

45-53 

F9.  2 

Target  altitude  hp  (m) 

54-62 

F9.  2 

Hufford's  integral  equation  (lath 
increment  As  (m) 

63-7  1 

F9.  2 

Measured  TDl  (psee) 

72-80 

FO.  2 

Measured  11)2  (psee) 

18, 

2-80 

Next  Target  Cards.  Use  formats  of 

Card  17. 

LAST 

CAmi 

54-62 

Blank.  Terminates  the  Input  Card  Stream. 

4.  1.  :i  STATION  CAHDS 

('.lids  14,  1.1,  and  16  specify  the  positions  of  the  .Muster,  Slave  1,  and  Slave  2, 
respectively.  Cards  1.5  and  16  also  contain  the  emission  delays  for  the  Slaves. 

4.1.4  T AIKiKT  C.MfIJS 


C.ird  17  contains  target  identifier,  position,  altitude  h^,  integration  path  in¬ 
crement  size  As,  and  measured  ITJl  and  Tl)2.  The  latter  pair  may  be  omitted. 


\n\  ■iiiiiiber  ')t'  :iiUliti(ini' 1  riMv  he  oiu'  l<>  ■.  ■."■1.  i  .u-d 

'.'.itii  Ji.  d  0  V. 'll  sioi)  I;!!'  Ill  iii-’ipi- '-ly  ifi'  p  li  jI' r  i ,  > 

M.  iik  bf  uiritM'lid  at'lr:'  :iip  lust  lafL’fl  laid. 

1.2  1  l.i|(r 

!'ii  ly 'la'  lltd'iM  rfads  ti>;  rain  a:  i  ii  i  nr  .la-  •  r  .j  a  I'i : : \i  Ir ,  .in.  i  ■  a  ijr.laiirp 

Diui.-iP  Ifni.n  Ibr  lalalKi.aP  taiir.  !  hr  ’.anr  rom.iin.-:  cinr  .‘ur  ii  dala  Irari  'nr  •-■vrr',' 
dO  .sr('  ')[’  Lalilndr  an  i  lon^iluur.  1- riah  triad  iS  packm  iiitn  -dtj  bii.s,  and  .siiirr  thr 
wiird  sizp  a\'ailahlr  on  a  C  lil'  GGOd  .  ; int [tiitrr  is  od  b:ts»  lU'O  data  triads  fit  into 
raali  GO -bit  word.  'I'hr  rie!:t  lia! ''-vvor.l  .ahvavs  rontaiiis  data  for  ihr  grtitri-aphic 
point  30-snr  north  of  thr  point  in  tim  Irl't  iialf-word,  as  indicatnd  in  labfn  2. 

Onn  database  rri  ord  rontains  sixlv  60 -bit  words,  for  a  total  o!  120  liata  iriads. 

.A  reronl  thus  spans  1"  of  latitiidr  from  south  to  north,  as  siiown  in  h  igurn  3(a). 

1  wo-hundrnd-fortv  surii  records  fill  a  database  zone  which  spans  2°  of  longitude 
from  west  to  east  in  30 -sec  steps. 

The  database  tape  must  contain  enough  zones  to  completely  cover  the  region 
.specified  via  the  map  corner  coordinates.  .An  e.'cample  of  the  ordering  of  zones  on 
the  database  tape  is  diagrammed  in  Figure  3(b).  Data  pertinent  to  zone  1  appears 
in  records  1  through  240,  to  zone  2  in  records  241  through  480,  and  .so  on. 


Table  2.  Database  Word  Packing 


Half 

Bit  Numbers 

Data  Description 

Least  Significant  Bit 

Left 

50-49 

Terrain  Elevation 

8  m 

48-41 

Impedance  Amplitude 

0.  001  n 

o 

o 

Impedance  Phase 

0.  00 1  rad 

Right 

29-19 

18-11 

/  Repeat  above  \ 

/  sequence  for  ] 

1  geographic  point! 

(as  above) 

10-0 

\  30  sec  north  / 

- 

1.3  I’rinted  Oiilpul 

T'hc  printed  output  consists  of  Normtil  Output  and  Error  Messages.  .Normal 
(jutput  includes  Control,  M:ip,  and  Station  parameters.  TDl  and  TD2  predictions, 
md  tile  End  AU'Ssage.  Figure  4  shows  a  sampb-  normal  output  listing.  Error 
■rio.ssiigcs  mav  appear  in  this  listing  in  the  event  of  eith-'r  fatal  or  nonfatal  errors. 
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RECORD  240 


TO  NEXT  ZONE 


WORD  ) 

POINT  120 

POINT  120  J 

POINT  120 

-  48 

59'  30" N 

60  ^ 

POINT  ng  ' 

POINT  119  , 

POINT  119 

-  48 

59'  00''N 

WORD  i 


FROM  PRIOR  ZONE  • 


RECORDS 
3  TO  239  j 


POINT  7 

1 

POINT  2 

POINT  1 

POINT  1 

48  00'  30"N 
48  00'  00"N 


Figrure  3(a).  Sample  Database  Zone  Structure,  Showing  Ordered  Sequence 
of  Database  Words  and  Corresponding  Geographic  Coordinates 


N.E.  MAP  CORNER 

54  00'  OO  'N 
—  53  59-  30"  N 

53  00'  OO-N 


—  48  59'30'N 


S.W.  MAP  CORNER  ■ 


Figure  3(b).  Sample  Map  Structure,  Showing  Ordered  Sequence  of  Zones  and 
Corresponding  Geographic  Coordinates,  The  coordinates  of  the  southwest 
and  northeast  map  corners  determine  the  number  of  zones  in  each  direction 
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0.00 
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4 
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4 

4b,3b4F 

)04ft.0C 

10,00 

0.00 

0.00 
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TAOPif  T 

e 
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4 
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TACGFT 

ft 
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3J 

•t 

4? 
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0,00 

10.00 
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13S01.J0 
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T 
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24,S44F 

0.00 

10.00 

0,00 

0.00 

»  3440.03 

?ft?53.69 

13444,83 

2b2S3.ft4 

T40GFT 

f 
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4 

S3 
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13484. OS 
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4 

44 
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134Hd,00 
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TADGf T 

10 

44 
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4 

? 
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0,00 

10.00 

0,00 

0.00 

13ftftrt,lO 
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13be8.l4 

2bSlb.lft 

Tf PMt^ATOP 

0 

(> 

o«ooo 

0 

0 

0«00y 
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0,00 
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ISPUT 

DATA. 

(b) 


Figure  4.  Sample  Normal  Output,  (a)  Control,  map,  and  station  parameters 
corresponding  to  input  sequence  01  Figure  2.  (b)  LORAN  time -difference  pre¬ 
dictions  for  target  cards  of  Figure  2 


4.  3,  1  CONTROL  PARAMETERS 

The  control  parameters  are  f,  t? ,  a,  a^,  b^,  c,  F^,  a,  electric/magnetic  field 
control  switch,  and  rough/smooth  earth  switch.  These  were  defined  in  Sections  2 
and  3.  Also  included  in  this  listing  are  the  calculated  quantities  k  (wavenumber), 
a/ (V  (effective  earth  radius),  and  the  spheroidal  flattening 


f  =  1 


o 


/b 

o 


and  the  first  spheroidal  eccentricity  squared 


1  -  a  “/b  ‘ 
o  o 
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4.  2  MAP  PARAMETERS 

The  Map  parameters  are  the  southwest  and  northeast  map  corners  defined  in 
Section  4.  1.  2. 

4.  3.  .3  STATIONS  PARAMETERS 

The  station  parameters  are  the  positions  and  emissions  delays  of  the  Master, 
Slave  1,  and  Slave  transmitters  defined  in  Section  5.  1.  3. 

4.3.4  T  D 1  A  N  D  T  D  2  P  HE  DlC  TIONS 

l  or  each  target  there  is  one  printed  line  which  includes  identifier,  latitude, 
longitude,  altitude,  numerical  integration  step  size,  and  measured  TDl  and  TD2 
defined  in  Section  5.  1.  4.  Also  included  are  the  computed  TDl  and  TD2  and  the 
time  delay  errors  (computed  minus  measured)  in  microseconds. 

4.3.5  END  MESSAGE 

The  message  "End  of  Input  Data"  indicates  normal  run  completion. 

4.3.6  ERROR  .MESSAGES 

.‘\  fatal  error  is  one  from  which  there  is  no  reasonable  recovery.  Hence,  the 
program  must  abort.  The  user  merely  corrects  the  fault  described  below  and  re¬ 
runs  the  program.  On  the  other  hand,  for  a  nonfatal  error  the  code  either  takes 
remedial  action  or  it  abandons  the  current  target  and  advances  to  the  next  one. 
Error  messages  are  grouped  below  under  the  routine  which  generates  them.  The 
message  number  refers  to  the  FORTRAN  statement  number  of  the  message. 

4.  3.6.  1  Program  HECTIC  Error  Messages 

Message  9900  (Fatal) 

".  .  .Station  lies  outside  map  area.  Run  aborts.  "  Fault:  Input  station  coor¬ 
dinates  exceeded  the  boundaries  specified  on  the  map  input  cords.  All  three 
transmitters  must  lie  within  the  map  region. 

Message  9910  (Fatal) 

"Error  in  data  for  .  .  .  coordinates.  Run  ends.  "  Fault:  Map  or  station  coor¬ 
dinates  on  input  card  were  not  correct.  Either  the  specified  latitude  degrees 
exceeded  90,  longitude  degrees  exceed  180,  minutes  exceeded  59,  seconds  exceeded 
59.999,  latitude  symbol  was  neither  N  nor  S,  or  longitude  symbol  was  neither 
E  nor  W, 

Message  9912  (Nonfatal) 

"Error  in  end-point  coordinate.  Calculations  for  this  path  have  been  deleted.  " 
Fault:  Target  coordinates  on  input  card  were  not  correct.  For  diagnosis,  see 
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message  9910  above.  Action;  The  program  advances  to  the  next  target. 


Message  9916  (Nonfatal) 

"Error  detected  during  conversion  of  coordinate  for  point  — from  radians  to 
(alphameric)  degrees-minutes-seconds.  "  Fault;  A  point  on  a  station-to-target 
path  was  not  correct.  Either  latitude  exceeded  324000.  0005  seconds,  or  longitude 
exceeded  643000,  0005  seconds.  Aclion;  The  program  advances  to  the  next  target. 

Message  9922  (Nonfatal) 

"Target  lies  outside  map  area.  Get  next  target."  Fault:  Input  target  coor¬ 
dinate  exceeded  the  boundaries  specified  on  the  map  input  cards.  Action:  The 
program  advances  to  the  next  target. 

4.  3.  6.  2  Subroutine  GEOPTS  Error  Messages 

Message  9800  (Fatal) 

"Subroutine  GEOPTS  called  with  end-point  latitude  out  of  acceptable  range 
(-7t/2,  ;r/2).  Latitude  of  A  =  ...  radians,  latitude  of  B  =  ...  radians.  Program 
execution  terminated.  "  Fault;  GEOPTS  encountered  a  bad  latitude. 

Message  9802  (Fatal) 

"Subroutine  GEOPTS  called  with  end-point  longitude  out  of  acceptable  range 
(-2n,  +  2n).  Longitude  of  A  =  ...  radians,  longitude  of  B  =  ...  radians.  Program 
execution  terminated.  "  Fault:  GEOPTS  encountered  a  bad  longitude. 

Message  9803  (Fatal) 

"Geodesic  path  includes  a  geographic  pole— subroutine  GEOPTS  cannot  handle 
this  case.  " 

4.  3.6.  3  Subroutine  GETELV  Error  Messages 

Message  9700  (Fatal) 

"Database  has  nnn  records.  Increase  array  MSINDX  and  variable  INDDIM  to 
accommodate.  Run  Aborts.  "  Fault:  Subroutine  GETELV  computed  the  number 
of  records  which  the  database  tape  must  contain  in  order  to  be  compatible  with  the 
specified  map  region.  This  number  exceeded  the  number  of  mass -storage  record 
pointer  locations  available.  Consequently,  the  database  could  not  be  loaded. 
Remedy:  the  user  should  increase  the  size  of  array  MSINDX  to  nnn  +  1  and  reset 
the  value  of  INDDIM  to  nnn  +  1.  Array  MSINDX  is  in  blank  common  and  INDDIM 
is  in  a  data  statement  in  routine  GETELV. 


Message  97  01  (Nonfat al) 

"Lat,  Ion  requested  are  outside  map  region  —  no  elevation  r-eturned.  "  l-'ault: 
The  program  autoniatically  iiie  rements  any  given  path  bj  two  stepsizes  bi.'yond  the 
tii'get  loeatioii.  Con.seqijentl,\'.  tlic  p.-.th  foi- a  target  that  initially  was  within  tiie 
map  region  ended  up  erossing  a  map  boundary.  Action:  The  program  proceeds, 
but  tile  results  tor  tliis  target  ar-e  questionable.  The  user  m:i>  try  a  smaller  step 
si/a  ior  this  target. 

Me s sage  “7Q2  (.\onfatal) 

"Index  requested  c.xceeds  scanlines  generated,  thus  previous  elevation  will  be 
used.  "  !•  ault;  The  record  pointer  returned  by  subroutine  INDEX  to  DETELV 
exceeded  the  number  of  records  on  the  database  tape.  Action;  The  program  pro¬ 
ceeds,  but  the  results  for  this  target  are  (Questionable. 

•!.:•!. 6. 4  Subroutine  CNEVKEN  Error  Messages 

Message  9602  (Eatal) 

"There  are  not  enough  points  in  the  given  array.  ”  Fault;  CNEVKEN  tried  to 
interpolate  a  function  which  was  defined  over  an  insufficient  range  of  values. 

Message  9600  (Fatal) 

"The  X  values  are  not  arranged  in  ascending  order. 

I  =  .  .  .  ,  X(I)  =  .  . . ,  J  =  . ...  A(J)  =  . 

X  F(X) 


Fault:  CNEVKEN  interpolates  the  funxtion  F(X)  over  array  X  from  a  given  function 
F(A)  over  array  A.  Both  X  and  A  must  be  in  ascending  order.  The  X  array  was 
not  in  order. 

4.  3.  6.  5  Subroutine  INDF  Error  Messages 
.Message  9500  (Fatal) 

"In  INDF,  distance  is  zero  or  negative,  XX  =  ..."  Fault:  INDF  tried  to 
calculate  the  induction  field  for  a  meaningless  distance  stored  in  variable  XX. 
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r>.  PROGRAM  IJSTING 

•MODULE  1.  4/?2/eo 

•DECK  HECTIC 

program  hectic  (  lNPitT=128,OUTPUT  =  l?HtTAP(-  1.1=  i £-8,  TAPE ?=!?«) 

C  INSTRUCTIONS  FOR  USE! - 

C  ATTACH.IN, INPUTCARDS, IO=TICHOtMR=l . 

C  ATTACHtHECTIC, ID=TICH0«MR=1. 

C  FTN,I=HECTIC,EL=E,OPT=1.L=0. 

C  VSN.TAPE13=CC3628. 

C  LABEL.TAPE13.R*NORlNr,,L  =  DATAHASE. 

C  LGO.IN. 

C  END  OF  INSTRUCTIONS. 

C  PROGRAM  DESCRIPTION! - 

C  MODULE  1  READS  INPUT  CARDS*  CALL  MODULES  2  AND  3*  AND  PRINTS 
C  FINAL  RESULTS. 

C  ROUTINES:  CORDMS.CORRAO.GEODI .HECTIC. 

C  MODULE  2  GETS  ALTITUDES  AND  IMPEDANCES  FOR  A  STaTION-TO-TARGET 
C  PATH,  READS  DATAEIASE  TAPE  IN  FIRST  PASS  ONLY. 

C  ROUTINES:  GEOPTSfPCOOPn.GEORET.OETELV, Index, INT, UNPACK . 

C  MODULE  3  COMPUTES  THE  SECONDARY  PHASE  FACTOR. 

C  ROUTINES!  GANG, CNFVKEN,6E0M, GROUND. FLEAF,INDF*lNEfci2E* 

C  OMCOS.WERF . 

C  files  USED! - 

C  INPUT  -  FOR  INPUT  LaRDS. 

C  OUTPUT  -  FOR  RUN  STATUS,  ERROR  MESSAGES  AND  DERUGGING. 

C  TAPE2  -  FOR  MASS  STORAGE. 

C  TAPE13  -  FOR  DATABASE  INPUT  TARE, 

C  REVISED  BY! - 

C  ELI  TICHOVOLSKY,  ARCON  CORP..  WALTHAM,  MA.,  4/PS/80, 

C  ORIGINAL  BY! - 

C  J.  R.  JOHLER.  ET.  AL. 

C  ENTRY  POINTS! - 

C  PCOORD  IS  IN  GEOPTS,  GEOM  IS  IN  GROUND. 

C - 

DIMENSION  TPW(3) *EO(2) .TO (2) . lOENTF (4 ) , TDME Ab ( 2) , TDERP ( 2) .TOA (3) 
DIMENSION  LDUM(3.4) ,A0UM(3,4) •RLA(3) ,RLO(3) ,L AT A ( 3) ,LON A ( 1) 

COMMON  /Z0TA/ARRAY( 15) 

COMMON  /PIS/TWOPI.PI.HAFPI.ORTPl 

COMMON  /MAP/LATSW,LATNE*LONSW,LONNE 

COMMON  /path/rlata.rlona.rlatb.rlonr.raza.razb.sbm 

COMMON  /CITCT/ITCT,LTOP.NPTCT,ZXMTR.ZMIN,AlTsw(3) *T5W(3) 

common  LAT(999) , LON (999) .MSINOX (5000) 

equivalence  (ARRAY(2) .ETA) , (ARRAY ( 1 1 ) .FLAT) 

equivalence  (ARRAY (12) .ALPHA) . (ARRAY (14) .EOHH) , (ARRAY (15) .SMOOTH) 
equivalence  (APRAY(I) .FKHZ) . (ARRAY! 13) .RAD) * (ARRAY (3) .SPHpaD) 
equivalence  (ARRAY(4).A0)*(ARRAY(5),B0) ,(ARRAy(9) .C) 
equivalence  (ARRAY (6) ,FL) . (ARRAY (7) .ESO) . (ARRAY (10) .WAVE) 

C  SEmI-MAJOP  and  semi-minor  axes  (METER)  OF  THE  INTERNATIONAL  SPhEROIO. 
DATA  A0,B0,SPHRAO/63783e8.O»  6356911.9461,  6367390./ 

DATA  LINE, ITCT, PI/50,  1.  3.1415926535690/ 

DATA  FKHZ »ETA,C/1 00., 1.000338,. 2997925/ 

DATA  flat , ALPHA ,EORH, SMOOTH/ 1OOO.,.05,1.,O./ 

DATA  WEST,SOUTH,LWEST,LSOUTH/lHW,lHS,lHW,lHb/ 

TW0P1=2.*PI  $HAFPI=PI/2.  $ORTPI*Pl/4. 


READ  PPOCES'^ING  CONTROLS. 

1=0 

20  READ  9919,I,AVALUE  5IF  (I  .LE.  0)  GO  TO  21 


22 


apwAY  (  n  =6v''H'^  >'■<  Tfi  ?<) 

^■1  CdMTIMiK 

c  pan=e  FFI-CT  I VF  FAPTi-  paplus.  wav*- =.;AVt''U)^'i'F  P  1Fi  «iP  at  F  AWTH  SMwFaCf  . 

PAD  =  SPHkAL)/Al  PHA  $wavF  =  T -/OP  I  *f  |\  H7  T  a  *  1  .  t -6/C 
C  FL  =  SPHEROir)AL  FLATTFiiI^G.  t‘^n=FlQST  htCF  f.'TP  I C  I  T '  S^jUApFO. 

FL  =  1 .  -  HO/an  its'i  =  1.  -  (MO/AO)«*r 
PRINT  R9?0.  ( ARRAY (  I )  » I  =  1 .6) 

PRINT  RS21,  ( ARRAY ( I ) . I =7«1S) 

PRINT  <5'i25 
C  READ  RAP  CORNERS. 

00  911  1=1,2 

READ  RROr,  lOEF  TE  (  I  )  .LOUM  (1,1)  ,LL»UR  (1.2)  .ADuE  (  f  ,  1  )  . .'  f):  .R  (  I  .  P )  . 

1  LOURd.O)  .LOURd.A)  .A0'IM(1,3)  .ADL'E  d.A) 

PRINT  PROP,  I  OF  NTF  (I  )  ,  LOOM  ( I  ,  1  ) ',  LOOM  d  .  2  )  .  A(.UEi  (I  ,  1  )  .  AOOR  (1  ,?)  , 

I  LOUR (1 .3) .loom (I ,4) ,ADUM( I ,3) « ADUM (1,4) 

MIS  =  0 

call  (  ORRAI)  (Rl.  A  (  1  )  .LOUMd  ,  1  )  .LOOMd  .  2  )  ,  AOOM  (1  ,  1  )  ,  AUUMd  .2)  .O.^l  S) 
call  CORRAO (Rl  0(1) .LOOM ( I ,3) ,L!)UM(I ,a) , AOUM  d ,3) , AOUM ( I ,4) , 1 ,M1S) 
IF  (MIS  .FO.  0)  on  TO  911 

c  IF  MAP  Coordinate?  are  rad,  ahort. 

PRINT  9910,  IDFNTFd)  SSTOP 
911  CONTINUE 

LATSw=LnuM( 1  , 1  ) *10000*LOUM(l ,2) *100*IF I  X ( ADOMd ,1 ) ) 

IF  (AD(jMd,21  .EO.  SOUTH)  LATSw=-LATSL 
LONSW=LnuM(l ,3) *10nOO*LDUM(l ,4) *100*IF I  X (ADUE ( 1 ,3) ) 

IF  (A()l/Md,4)  .Elj.  RFST)  LONSws-LONSR 

LATNE=LOUM (2,  1  ) *  1 0000 *100” ( ?, 2) *  1 00* I F 1 X ( ADLN ( 2, 1 ) ) 

IF  (A0UM(?,2)  .EO.  SOUTH)  LATNF=-LATNf 
L0NNE=L0uM(2,3) *10000*LOUM(2.4) *1 00* if  1 X ( AOCE  (2,3)) 

IF  (ADUM(2f4)  .Eu.  WEST)  LONmE=-LONNF 
C  READ  transmitter  COORDINATES  For  MASTER,  ?L AVE  1  ♦  ANO  SLAVF2. 

PRINT  9926 
00  9S0  1=1,3 

READ  9909,  IDE  NTF  (1 )  ,  L  OUN  (  I  ,  1  )  .LI)UM(  1.2)  .uDl.N  (1,1)  ,Al)UMd  ,2)  , 

1  Lt)UNd,3)  ,LnuX(I,4)  .ai  U-M  1 ,3)  .af)llM(I,4)  .datum 

print  9909. 1  of  f  TF  d  )  ,LDU-»(  I  .  1  )  .LDUM  (1,2)  .aDDM(I  ,  1  )  .  ADUN  (1.2), 

1  LUl'M (1.3) .LDUM( I ,4) . ADUM  d .3)  , ADUN (1,4)  .DATUM 

IF  d  .YE.  2)  Fi)(I-l  )  =  datum 

MIC.  r  0 

call  C(  RRA('  (R(  O  ( I  )  .LOUMd  .3)  .loom  d  .4)  .  Al'UM  d  .3)  ,  ADUM  (I,4),1,MIS) 
CALL  CORRAO  (RL  A  (1  )  ,  LOIlM  ( 1  .  1  )  .LDUM  (1,2)  ,  AIJUM  d  ,  1  )  ,  ADUM  (I  ,  2 )  .  0  ,  m  I  s  ) 
IF  (MTS  .EO.  0)  00  TO  940 
C  IF  STATIQf  COOPDIF'ATFS  ARE  PAD,  AHURT. 

PRINT  991  0.  IDFliTF  (I  )  SSTOR 

940  LATA (I ) =LDUM(I ,  1  ) *  1 0000 .Ld  im (I , 2 ) « 1 PO ♦ 1 F I  X (“RUM (I , 1  )  ) 

IF  (ADUMd,2)  .EU.  SOUTH)  L  AT  A  (I  )  =-l.  A  T  A  (I  ) 

LON  A  (  I  )  =lDIIM(  T  ,  3)  »1  0000  .LOllMI  I  ,  c )  «  1  00  ♦  I  F  I  X  (  ADUM  (1,3)) 

IF  (A0UM(I,4)  .EiJ.  /(EST)  LONA(n=-LnMA(I) 

IF  (LATA  (1  )  .CtE  .LATSi-i  .ANLi.  L  A  T  A  (I  )  .  LT  .  L  A  TF-E  ,AND. 

1  LONA (I ) .OF.LONSw  .and.  LONA(l) .LT.LOnNE)  GO  TO  9bO 
C  IF  STATION  LIFS  OUTSIDF  MAP  AREA,  ABORT. 

PRINT  990C, IDF NTF (I )  SSTOP 
9S0  CONTIMiF 

C-  -  -  -  - - - - - - - - - - - ----- 

c  start  LOOF'  on  targfts. 
s  continue: 

IF  (LINE  .1 T.  SO)  00  TO  SS 
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-d-fif.  TA>.(,fT  CfHiPU  U'uT  F /“U;  =  A  I  wHO'dHf  -  T  T  ^LTITIIlOK  (M)  AMOVfc. 
(,  MKAti  St  A  LFVKL  if  M(>N-;fpi).  All^l.S  =  PATH  INCPFPtNT  SI7F  (Kw,). 
hvtfl'  (H)  .LAT'.ti)H.LArPirjH.SFCLATH*LATlDF'. 

.  I  ONPF  GP  .  L  ON'-'  I  >'!■■<  »  SE  CLONB  ,LOf.'  I  Oh  .  7  P  1  ,  ADE  L  S  .  T  DmE  As 

PPIST  wSOg,  IDFEiIE  (A)  ,l  ATrjFGri,LATMlNH.SECLATP*LATII)t-i, 

1  LOMPE(,H,LO^JMIN^-.SFCL<l^lH,LO^•li'H.7Ml^J,ADELStTI)Mt  AS 

LINE  =  i,  I\F  *  1 

L  STi)P  IF  F  AJH  PiCWEpF  mT  SI/f  is  7F  WO  OS’  LESS. 

I  F  (  AOt  IS  .(  E  .  0  .  )  TO  in 

Pis  =  u 

Call  C0PPA(.  (PI  ATh.L  AT’.jF  C'H  ,L  ATHIMH  ,SECL  ATH  .LA'lDHt  dtMl  5) 
call  f  pPPaO  (pi  OH'-  ,L  OfjliF  lip  .LO^^Ml^H  .SECLONH  ,LlA„  lOe,  1  ,MI  S) 

IF  (PIS  .Eg.  C)  (p.  rn  IS 

C  IE  TAPr,ET  COOSPlNATtS  APF  HAI.,  GET  ANOThFP  TAPGLT. 

PWlgT  SPlie  EL  l''  F  =  l  I'lF*"' 

(iO  TO  s 

f_  ___________________________  ______ 

L  STAPT  LOOP  ON  .1  S  T  A  T  I  OE'-TO- T  APC-F  T  PATHS. 

IS  continof 

uLATA=PLA (TTCT) 

PLONAiPLO ( I TCT I 
LTOP  =  (i 

ALTSW(ITCT)  =  TSw(ITCT)  =  SBKMS  =  0. 

C  SKIP  This  path  if  it  is  too  smopt. 

IF  ( (AESJPL ATA-plath)  .LE.  AOELS/SPhPAD)  .ANU, 

1  ( APS (PLONA-klONB)  .LE.  AOFLS/SPHPAD) )  GO  TO  16 

c  Compute  sphe.pciioal  baseline  and  fopwaro  amo  baca  a7Imuths, 
call  gfopts (n .n .wlatp*plonp»pa2P) 

C  COPPUTE  SPHE'woiOAt.  hASELINE  AEjO  FOPWAWO  AND  HACA  A?rMUTHS, 
call  GEOni (PLATA, RLONA,PA7AS.SBMS,PLATB*WLONB,RAIBS) 
SPkMS=SPMS*.001 
SFiKP  =  SHP«l.rE-3 
OSkm=nptCT=SHkp/ADELS*1 .0 
osp  =  sFiM/nsKP 
OSkp  =  OSP*1.0E-i 

C  SET  tramspittep  cooroimates  as  point  number  1. 

LAT(1)=LATA(ITCT1  SLONd  )  =L0NA(ITCT) 

C  ADD  2  POINTS  BEYONO  RECEIVER. 

NPP^=NPTCT*2 
SPM  =  0.0 

C  compute  latitudes  AMU  LONGITUIYES  ALONG  HuFFOPD»i>  BASELINE. 

00  100  IP=1.NPP2 
SRm  =  SPM.TiSM 

CAl  L  PCOOPO (SPM, PLATO, RLONP.RAZP) 

MIS  =  0 

call  CORDMS (RLATP,LATDEGP,LATMINP,SECLATP,L*'T I  dp, O. MIS) 

call  COPDMS (PLONP,lONDEGR,LONMINP,SECLONP,LONIPP, 1 ,MIS) 

IF  (Mis  .F(J.  0)  GO  TO  3S 

C  IF  PATH  COORUINATF  IS  RAO,  GET  ANOTHER  TARGET. 

PRINT  9916. IP  slINE=LINE*4 
GO  TO  S 

35  COMTINmF 

LAT ( IP, n  =LAT0FGP*1 noOO*LATMINP*100» IF  IX (SECLATP) 

IF  (LATIDP  .EO.  LSOUTH)  LAT ( IP*1 ) =-LAT ( IP*1 ) 

LONI IP*1 ) =l  ONOFOPO] OOnO*LONMlNP*100,IFI X (SFCLONP) 

IF  (LOMOP  .EO.  LWEST)  LON(lP*l)=-LON(IP,n 
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100  continue 

C  100  1  POINT  FOP  TPANSMITTEP. 

1P=NPTCT»3 

C  IF  TARGET  I'",  OUTSIOE  MAP  AREA.  GET  ANOTHER  TARGET. 

IF  (LA',  (  IP)  .6E.LATSW  .AND.  L  AT  (  I P  )  .  LT  .L  ATNt  .AND. 

1  LON  (IP)  .GF.LONSv*  .AND.  LON  (  IP )  .LT  .LONGE  )  GO  TO  8^0 
PRINT  99?P, lOENTF (4)  $L I NE =L I NE *  1  tGO  TO  S 

C  COMPUTE  altitudes  AND  IMPEDANCES  ALONG  HljFFORO'b  BASELINE. 

0SO  call  GEORET (OSKH, IP) 

C  ADJUST  NUMBER  OF  POINTS  FOR  ROUTINE  INtOPE. 

NPTCT  =  NPTCT  ♦  1 
C  COMPUTE  HUFFORO'S  EQUATION. 

CALL  INEQ2E 
16  CONTINUE 

TPW(ITCT)  =  ETA*SBKMS/C 
ITCT=ITCT*1 

IF  (ITCT  .LT.  4)  GO  TO  IS 
C  END  LOOP  ON  3  STAj I ON-TO-TAR6ET  PATHS. 

C-  -  -  -  -  -  - - - - - - - - 

C  PRINT  FINAL  RESULTS. 

00  9  1=1,3 

IF  (LTOP.£<J.O)  TOA(I)  =  TPW(I)  ♦  TSWII) 

IE  (LTOP.EO.l)  TOA(I)  =  TPW(I)  ♦  ALTSW(I) 

9  continue 

DO  710  I  =1,2 

TD(I)  =  ED(I)  ♦  TOA(I*l)  -  TOA)l) 

710  TOEPRd)  =  TO(I)  -  TOMEAS(I) 

PRINT  4923,TD,TDERR 
C  GET  NEXT  TARGET, 

ITCT=1 
GO  TO  5 

C  END  LOOP  ON  targets. 

C-  -  --  --  --  --  --  --  --  --  --  --  --  --  --  --  --  --  - 

10  PRINT  9918 
STOP 

9900  FORMAT  (»  *A10*  STATION  LIES  OUTSIDE  CAP  AREA.  RUN  AkoRTS.o) 

9909  FORMAT (2X. A  10,2 ( IS, I 3.F7.3.A1 ) ,4F9.?) 

9910  FORMAT (*0  ERROR  IN  DATA  FOP  »A10*  COORDINATls.  RUN  ENOS.o) 

9912  FORMAT  ( //1 1 X,  anttoo,  9x,  29HERROR  IN  FND-POlNT  COOROI'-aTE,  1 0  < .  "(H 

♦*«o/l lx .SAP***  calculations  for  This  path  have  been  DELETf)'  o>«) 

9916  FORMAT(//2x,  T RROR  DETECTED  DUPING  CONVERSION  oF  CDORi'riAT 

♦  ES  FOR  POINT,  IS.  SH  ***/2x,  ?H««,  7X.  S2hE  kOm  RADIANS  Tr  (AU-,-Ai' 

♦  ERIC)  OEGREES-MINuTES-sECONDS,  7x,  3i-,o«*/?x,  3h**«,  h6X  , 

.  3H***,  7X,  S2HFUPTHEP  CALCULATIONS  TOP  THIS  PATH  HAVE  HETn  TEL^T*- 

*D.  7X,  1H*«*) 

9918  EOPMATC'O  end  OF  INPUT  OATA.o) 

9919  FOPMA I ( IS,E 1 n.O ) 

9920  FORMAT  (IHI//1S(2H  -)*r,nNTPOL  P  AS  A..t  Tf  PS  "  1  S  (  2H-  )/ 

AE  ?2 . 1  .  SX  «♦  ARPA  Y  I  1  )  =  FKh/  =  FWFQUETiCY  (KH^)-/ 

1E2,’.  iP  .SX«^  AppA  Y  (2)  =  FTA  =  AIR  INDEX  OE  PETpALTION  AT 
?E?2.  l4.Sx»,APPAY  (3)  =  SPHPAO  =  sphEPICAL  EAEIH  PAI'IUG  (iX)o/ 

4R22. 14,  Sx*>.  APPA  Y  (4  )  =  AO  =  spheroidal  EApTH  tr;uATUPIA|.  RAuiiic;  (  x- )  e 
S/E22. 14  ,SX«-ARRAY  (S)  =  HO  =  SUHEpOIOAL  E.kTh  polAL  PAl'ilUS  (M)o/ 
7>E??.14,SX<*  ARRAY(H)  ::  F  L  =  SPHfPOIOAt.  ELATTL.ING*) 

9921  format  (E.2?.14,jX0  array  (7)  =  E  Su  =  SPHE  PO  K)  A|.  FIRST  EtL^^TPlLlTY 

2  SUUAPrn**/ 

2E2?.14,SX*  ARPAY(8)  = 


NOT  USED*/ 


■*«!>- 


H  .SX**ARPAY  (<?)  =  vacuum  SPFtU  OF  L  U)HT  (KM/USfcC)*/ 

5E22,l4,5X*  ARRAY(IO)  =  WAVf  =  IvAVEfYl/MFe  R  IN  AIR  AI  EARTH  qilR 

7EACe*/ 

6e??.  14,5x»*array ( 11 )  =  Elat  =  flat  farth  factor  (H)*/ 

BE22, lA .bX** array ( 1? )  =  Alpha  =  VERTICAL  LAPbt  FACTOR*/ 

9E2?.1a,5X*  APRAY(13)  =  EEFECTIvf  EARTH  RAOIUb  (M)*/ 
1E22.1a,5x»*ARPAY( lA)  =  EORH  =  o.  FOR  H-FIFLR*  NON-O.  FOR  E-FIELO*/ 
9E22.1a,SX**ARPAY(1S)  =  SMOOTH  =  1.  FOR  SMOOIh  FARTH,  NON-1,  FOP  RO 
2IIGH  EARTH*/ 

327x**  ITEMS  mapkeo  with  ♦  may  he  chanoed  via  Input  cards.*/) 

9922  format  (*  *Alfi*  TARGET  LIES  OUTSIDE  MAP  APtA.  GET  NEXT  TARGET.*) 

9923  format  ( 1h. ,7qx .aF9.?) 

992a  format  (1H1//??(2H  -)*lORAM  TImE  OIFFEREnCF  HR£ 01 CT I0NS*2 1 ( 2H-  )/ 

12X*IOFNT1FIER«7X»LATITUDE*7XoloNG1TIIDE  altitude  STEPSIZE  _ MEASUP 

2En(USFC)-  — COMPUTED (USE C) -  ERROR (DStC ) --“  «/ 

312X.2(3X»  (DEG-mIN-sEC  )“)./X*(MF:TER)*SX*(KM)*J(6X*T()l«Gx*Tn?*)/) 
9925  format  (//8(?h  -)»^.ap  p  ARAMFTFRS*”  (  2H-  )/ 

12X*IDENTIFIER*7X*L ATlTUDE*7X»LONGlTUDF*/ 
212X.2(3X*(nEG-MlN-SEC)*)/) 

9928  FORMAT  (//11(?H  -)*STATI0F'  P ARAMf, TEPS*  1 1  ( 2H-  )/ 

1?X*IDENT1E1ER*7X*LATITUDE*7x*LONGI7UDE  EMISSION  DELAY*/ 
212X,2(3X*(DEG-mIn-SEC)») . 3x* ( usEC ) */) 

END 


•deck  coprao 

SUBROUTINE  COPRAO (RCOR. IDEG. IMJN.SEC. ID. IS.  ItRR) 

DIMENSION  lOS(A) 

data  (IDS=1HM.1HE.1HS.1HW) 

ISS  =  IS 

IF  (ISS)  10.5.15 
5  IE  (in.EO.IOSd  )  )  GO  TO  25 
IE  (ID.ED.IOSOI  )  GO  TO  30 
10  IFRR  =  1 
RETURN 

IS  IF  (Iss-1)  20.20. in 
20  IF  ( ID.EO. IDS (2) )  GO  TO  25 
IF  (IO.EQ.IDS(A) )  GO  TO  TO 
lERR  =  1 
RE  TURN 

25  SIGN  =  1.0 
GO  TO  35 
30  SIGN  =  -1.0 
35  IF  (IDEG-IaO)  AO. AO.  10 
AO  IF  (Imin-ho)  AF.m.lO 
A5  IF  (SEC-60.0)  CO. 10.10 

50  RCOR  =  SIGN*  (FLOAT  (  IDEG)®  ( I  .7A532925199A33E -<;)  ♦FLOAT  (  IMIN)« 
♦  (2.90HS82CH665722F,-A)  ♦  St C*A . BaS  1  368 1 1 09536^-6 ) 

RETURN 

END 
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DECi^  COWUMS 

SUHBOUT  IMf-  CORnMS(PCOB.  lUtG.  IMIN,►^EC.  It  .IS.  IEBM) 

DIMENSION  IDS (4) 

data  (IOS=lHN.lHE.lHS,lhV ) 

PAkO  =  PCOP 

SEC  =  Arts  (RANG)  «-20f^a64.80b?470'i6 
ISS  =  IS 

IE  (ISS)  S, 10,15 
5  lERR  =  1 
RETURN 

10  IE  (SEC-324000.0005)  25.5,5 
IS  IE  (ISS-1)  20.?0,S 
20  If  (SEC-84«000.000S)  25,5.5 
25  IE  (PANG)  30,35,35 
30  ISI  =  2 
GO  TO  40 
35  ISI  =  0 

40  IDEG  =  SEC/3600.0 

IMIN  =  SEC/60. 0-60. 0"EL0AT (lOEG) 

SEC  =  SEC-3600. U*fLOAT ( IDEG)-60.0*ELOAT ( IMIN) 
ISEC  =  SEC 

SEC  =  SEC-ELOAT (  ISEC) 

IE  (SfC-O.RRSS)  60,45,45 
4S  SEC  =  0.0 

ISEC  =  ISEC*1 
IE  (ISEC-60)  60,50,50 
50  ISEC  =  0 

IMIN  =  IMIN*1 
IF  (IMIN-60)  60,55,55 
S5  IMIN  =  0 

I  deg  =  IOEG*l 
60  lCx  =  ISIMS5M 

riEC=EL(!AT  (ISEC)  ♦SEC 
ID=IDS(LDX) 

RE  TURN 
END 


SUt-'-dl'T  iNf  Gt'TUC'l.  6Ll»  St  ,  4l?.  A/t-l) 

C  Pi'wmsf  : 

I  CfiLCULATf';  tOMPiHfiTKiN  rOpw  ( S'H'Afi'.' t  lSfc,S). 

L  tOiiri-fPN  L  AT  I  TCJOFS  AfiF  WESTfSG  L  I  Tl/f/f  S  Awf  ’jEGATJVh. 

c  If  puts: 

C  =  (tOlHTIC  LATIT,M>^  (PfiOlAU)  ,,E  FWST  POINT. 

C  A(.  1  =  oFoni-'TK  LOnt-nimF  (sAI'M'-O  of  fjust  pojtt. 

I,  =  OEOOTTIC  LATITunF  (PftDlAD)  iiF  StCONt  HijlUT. 

L  AL^*  =  OtOf'ETU  LOTGITunt  (PADIAN)  Of  ST  C  (U  O  POJtiT. 

C  outputs: 

L  A71?  =  FOPWAPO  A7I’‘’uTh  (PAOIANi)  F pOm  FIPsT  pCjINT. 

C  A/<?1  =  HACK  A7IM.iTrt  (PAOIAT.)  FpOP  St  COS'I  PoiNT, 

C  S  =  otOOfTIC  niSTANCE  (T*t:TtP)  pFT/jEFM  P(  IHTS. 

COf'MOF  7PIS/TwnPI  .PI  .hAFPl  tOPTPI 
COMMON  /70TA7APkAY ( IS) 

Ou=APPAY(SI  tFL=APhAY(6) 

IF  (Af»S(Hl)  .r-T.  OPTPI)  GO  TO  1 
THETl  =  tan (Hi  1  *  ( 1 •  -  f  t  ) 

HfTl  =  ATANdHtTl) 

GO  TO  ? 

1  CHETl  =  1.0  7  TAN(BI)  /  (1.  -  FL) 

MtTl  =  ATANd.  /  CT-ETl) 

?  COMTINUE 

ALLl  =  AL?  -  ALl 
IF  (AL2  -  ALl  .EO.  0.)  A.  9 
H  ALL^;  =  0. 

GO  TO  3 
9  CONTINUE 

ALL?  =  AL?  -  AL  1  -  SlftN  (TSOPI.  Al  2  -  ALl) 

3  I‘  (APS(ALLI)  .GT.  APS(ALL?))  S.  h 

S  all  =  all? 

GO  TO  7 
b  all  =  ALLl 
7  CONTINUE 

12  IF  (AF.s(ALL)  .FO.U.  .OP.  ArS  (  ALL)  .EO.PI  .  OP  .  ApS  (  AL  L  )  .  F  '  •  T  -MP  1)  1  0  .  ]  1 

10  all  =  AhS(ALL) 

11  continue 

IF  (ApS(P2)  .GT.  OPTPI)  GO  TO  lb 
TPFT2  =  TAN(h2)  *  (I.  -  FI ) 

HFT2  =  ATANITPfT?) 

GO  TO  1 7 

16  CpET?  =  1.0  7  TAn(H2)  /  ().  _  FL) 

HtT2  =  ATANd.  /  CHFT2) 

17  Continue 

T 1  =  cos  (HFT  1 ) 

ScETl  =  SIN(pFTI) 

COT  T?  =  cos (DtT?) 

SpfT?  =  siN(eFT2) 

A  =  SrETl  'T  SPFT2 
H  =  CHFTl  *  CBFT2 
AP  =  shETI  »  CSfT? 

Mfi  r  CSFT2  ®  CPfTl 

cost.  =  COS  (  All  ) 

SINL  =  SIN(AlL) 

'FF-I  r  A  «  H  »  CnSl 

SPHi  =  SOFT  ((SINL  "  CPtr/)*®^  ♦  (iA  -  A-  U  cost)**?) 


2!1 


If  (bihi 

.f-f  . 

0.  .AMt'.  (.Pr’1  .(it.  I.l 

20.  2) 

WHl  =  ;.SI 

L  (  SPp 

1 ) 

If  (Fir.! 

.(-.T  . 

rPHI)  Uril  =  ACObiCPHl) 

oO  TCI  70 

?l 

If  (SPMt 

.Of  . 

0,  .A-iTi.  cpi-ij  .lF.  0.) 

20  .  P  * 

?? 

PHI  =  PI 

-  AblP  (OPrU 

If  (SPHi 

.OT. 

APf  (  CPfi  I  )  )  '■'H  1  =  ALC)S  (  r 

,PH  1  ) 

.•,f,  TO  30 

?.< 

If  (FPr-1 

•  LT. 

0.  .Asj,>,  C.PhI  ,1.  T.  IJ.I 

2a,  2b 

Phi  =  PI 

-  Abl 

f-  (SPhI  ) 

if  (Ai-b(SPPl) 

.l^T.  A4S(CP^^I))  ^nl  =  ) 

00  TO  3f' 

If  (SP--I 

.{  T. 

Cl.  .AN.J.  CPHI  .liP.  (/.) 

ir't' 

rJi-l  =  T41OPI  4 

1  SILiSP-'l  ) 

[f  (AOSibPHlJ 

.OT.  CPHI)  PHI  =  T^f'pf 

-  I  ) 

(III  TO  3  0 

^7 

CALL  FxIT 

JC 

COCiTlMl^f 

C  =  »  SiNL  /  S'-'HI 

f- 1  L  ®  fL 

CfiM  =  Ft  *  f  I? 

COKi?  =  n.S  »  K  ? 

cr'\'3  ::  SPHI  »  f.t-'r'I 
CDf.t.  =  Pt-I<**2  / 

CO’vS  =  CFIKU  *  C'-'r'I 

f  ^^  =  1  .  -  C  *  *  2 

l«ATlf)l  =  <!•  ♦  CjNl)  »  ^Hl  4  A  e  (CO'-l  "  S^HI  -  f  l.'N2  *  C0^!4) 

j  ,  t '•'  »  (-O.S  *  C  O^•l  *  (P-Hl  4  ';0N3)  4  (  «  CG^  S)  -  A  o  A  “  CCA!? 

2  »  C0^.3  4  ty  *  Em  4-  C0n2  «•  (0412^  *  (Gi-I  ♦  COW.!)  -  C  * 

3  00^3  <t  CP^'I**2)  4  A  4  f*'  *  C'1.%2  “  4  C'>M3  • 

s  =  kJATIOl  *  f-O 

IF  (S  .f  .  1  .t'-'*)  F  =  <■'. 

kATlOf'  =  Cr,M  *  Fi-I  -  A  *  C:')''J2  *  (S<'’4’I  4  A.  '•>  c.li"**) 

1  4  o.‘^  »  (■  M  «  ('.nN2  “  t-v.f'  <■•  4  CA'NO  4  A.o  *■  roAS) 

alaM  =  WAT  1 02  *  C  4  ALL 
SAl  AM  =  SlN(ALA") 

CALAM  =  CnS(Al,AM) 

CTA212  =  Wu  -  C.ALAK'  4  Ah 

CTA221  =  rA  4  CALAM  -  uA 

IF  (ALl  -  Al2  .tj.  0.)  3S, 

3S  A712  =  Cl, 

A/f-l  =  Cl. 

•j(i  TO  3 A 

3G  CTAZ17  =  rTA712  /  (SALAM  »  r-^eT-^l 
If  ( C T  A  7  1  2  .to.  A  .  )  S<4 ,  bF 
b<4  A712  =  HAfhl 

(j(j  TO  bh 
Sb  CONTIfiit 

A712  =  ATAt  (1.  /  rTA712) 

bh  COf.Tlr.iiF 

CTttZ21  =  CTA/21  /  (SALAM  *  c-hTl) 

IF  (CTA721  .tn.  .1.)  S7.  b« 
b7  A721  =  HAFPI 
00  TO  3A 
bH  C(K-TIfA)t 

A721  =  aTAM(1,  /  fTA7?t) 

3A  CONTIM'F 
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If  (Ail  .r.f .  o.  .Af'r>,  cta/j?  .ge.  o.) 

40. 

40 

A/ 12  =  a.M? 

60  TO  bO 

4l 

If  (A(.l,  .Of.  0.  .AMT.  CTA71?  ,LT.  0.) 

42, 

43 

4? 

A712  =  PI  *  A/1? 

bO  TO  bO 

43 

If  (All  .LT.  0,  .AfiD.  CTA712  .GE.  0.) 

44, 

45 

44 

A712  =  PI  ♦  A71? 

fiO  TO  bu 

4S 

A/12  =  T.IOPI  ♦  A/1? 

bO 

If  (AH  .Of.  0.  .Af'O.  CTu??!  .GE.  0.) 

46  . 

4  7 

46 

A7?1  =  pi  *  A7?l 

GO  TO  bl 

47 

If  (ALL  .OE.  0.  .AfO.  CTa/?1  .LT.  O.t 

4b, 

49 

46 

A7?l  =  T^OPI  ♦  A7?1 

GO  TO  bl 

44 

If  (ALL  .LT.  0.  .Aa'O.  CTA??1  .GE.  0.) 

b2. 

53 

A7?l  =  A7?1 

GO  TO  51 

b3 

A721  =  PI  ♦  A7?1 

bl 

CONTI  fJL'E 

=  AMOO  (A71i?«  TVOPI) 

IF  (A712  .LT.  0.)  A712  =  A712  ♦  TWOPI 

A/21  =  AMnO(A721«  TWOPD 

If  (A7?l  ,LT.  0.)  A7?l  =  A721  ♦  TWOPI 

Pf TUPM 

£F.D 


r 
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•MODULE  ? 

•DECK  UtORfcT 

SUPHOUTINt  6E0RET (DSKM.NP) 

COMMON  LAT(999)  ,L0N(999)  . 1  Ml h  ( BfiO 0 ) 

COMMON  /GRN01/7PP./0(ft) ,7(1009) ,7Pi)(0) .7P(1009) .XUlb) ,« (lno9) 

COMMON  /TE/ZPPl  (  1009)  .►'PT'; 

common  /SOROI/S  (999)  ,!)P(999)  ,01  (999)  ,LIM 

COMMON  /C1TCT/ITCT,LT0P,NPTCT,7*MTP,7min,alTS»  (3)  ,TSm:1) 

NPT5=NP 
S(1 )=0. 

X  ( n  =0. 

DO  1  I=1,NP 

call  (jETELV  (LAT  ( I  )  ,LON(n  .AAmP.EA/./ (  I  )  ) 

lE(I.GT.l)  S(I)=S(1-1)*0SKM 
OP(I)  =  AAMp  *  C0S(FA7) 

01(1)  =  AAMP  *  9IN(FA7) 

1  CONTINUE 
LLM=NP 

C  COMPUTE  i3ASEL1NE  TtkPAiN  ELEVATIONS  PeLATIVE  T(i  TWAK  Sl-'l  T  TEP  AITiTUOP. 
7XMTR  =  7(  1  ) 

Z(1)=0. 

00  2  1=2, NP 

c  FOP  no9,  convept  oistancf  to  MFTEPS. 

S (I )=S(  1) *1000. 

X(I)=S(1) 

7(n=?(I)-?XMTP 

2  CONTINUE 

C  COMPUTE  elevation  DERIVATIVES. 

NPM1=NP-1 
00  3  L=2.NPMl 
CALL  1NT(L»L> 

ZPPl (LI=7PP 

3  CONTINUE 
2P(I)=0. 

ZPPI (1)=0. 

2P(NP)=7P(NP-1) 

7PP1 (NP)=7PP1 (NP-1 ) 

PFTUPN 


SV't^'^ODT  l\f  lilt  f'P  t  ^  •'P) 

C(iv^'^)^  /P  K/T-()P  I  ,‘- !  t'-^M'l .  .-I^l 

iji  /PAJH/Pi  HI  f.  al  r- ,i.L‘ r,-.. -'t. 'i  . -u/ 

(Xv^v.f  //!■!  A/£Ut^iV  (  1‘-, ) 

Af  -  >■  S'-.  =  AwPAY(/)  i-CtS,.  -  \,  -  PS''. 

[^  (  A'.P  (-vt  f.T  A  )  -►■Si' wl  )  S.lf,.1(t 

s  IF  (  A‘'S  (  'L  AT.l) -p'AF  PI  1  IS.I.'.n 

l(j  Hwt^'l  ■*--i)n«PL''la*PLATH 

a>~ut.  F  1  iMi' ,  fiQP  SI.MP(  111  I  St  .'-.tOPTS  CALlF  ’  ■-•  IT'"  f  '  ' '.-Pf' I  J  t  i^aTIT.il 

♦  t.  <  UT  F- FA--L  F  PAS.jt-  ( -P' ) /i*.  .'^I  //F  ) // 1  3*  .  1  SP-I  M  r  -'•.  r.r  A  *, 

♦FP»y'*  pY-  SAti!A’',s  l.flli’.'  F  OP  't  =»  t  "  ,  ?  *  an  SA.ili'Y) 

F  •'.'  1  ^  T  p  S  u  1 

VfOl  FOPmAIi'-Ja,  pa^-P-f.r-.-Atp  t  >F  Csl  li'S  Tp  ,*.-lt.,ATF  Dl 
STOP 

IS  If  (  A‘- s  C'flc.Aa  1  -  f /'TP  1  )  P<'.p<'.PS 

f(i  If  (  AOS  (Y  l'T--'  ) -T  »r^' I  »  SO.jOttFO 

PS  fPtST  SHnp,Ptr.,A,PL  )A'- 

'toOP  fOppA7(lP0*  S‘JSPO''Tl‘'t  (tfOPTS  (IAllTi'  fITp  ffv)— POIjT  LONOITI 

♦  Ilf  OUT  of  ACCfPTAF.i  C  1.  ptvot  (-XPI  tPPl  )// 1G>.  1  THLO'-lFlTliOf;  of  a  *« 

♦  fs.p,  TIP  PAlJA'tS  Ll'f^JTO'f  fTf  r  fs*p«  pp*  PAi'IASS) 

PPIST  PFini 
S  T /p 

30  alonp=pu'^oa 
t-L')NP  =  PL.T'Fi 
If  (Ai  f  .jpFpi )  Is,  Ts,<.n 
3S  ALnNPiALOsp^T.f'Pi 
sn  IF  (PLOf'-*PI  )  AS.ss,'^f) 

SF.  HLf'f.P=PLoA!P*Tl-.OPl 
sn  SL  Ofj  =  -ILI'P.P'AI  'INP 

If  (PI  ('iY*Pl)  SS»70»7'y 
■'S  If  (eL‘''’'p)  t.3.sS.PF> 

M'  SLOf.FF  =  sLONw*T*OPl 
on  TO  PS 

SS  ALONp  ='  Al(.vNP-T«OPI 
GO  TO  S5 
70  PPIMT  vpO.3 

i,^03  FOpMATilHO,  (?PH  OFO.JfSlC  PATH  INCLUDES  A  GeoGSAF^K  polf  -  s  HP 

♦  OUTlhf  IjFoPTS  CAAif.OT  haOIaLF  THIS  CASf  ) 

PPINT  SHOl 

STOP 

7b  If  (SLC'm-PI)  pp.fo.eo 
SO  If  (SloM^-PII  SO.Pfi.SS 
8S  t<l  OAiP  =  SLONF-TS'uPI 
00  TO  P'j 

PC  Al  OOP  r  ALpTJP'TF'UP  I 
PS  =  PLAlA*PLATh 

PtTA  =  SIM0.p*S-H) 
pfTA  -  ( I  .n-fc  SO'Fpf  T  AoptTo  ) /Cf  SO 

mLAT  =  AO/ (MtT  A*S<  ;PT  (pf  TAsCFS'))  ) 

♦LO'.  =  SLOT. P-Al  OOP 
ALOf'  =  “iF  TA»ijLr’v 
0  =  SIf’(A|,iiM 
pLAT  =  ''oSiPLATH 
ALAT  =  COS(PlATA) 

SA/  A  =  O-Itbl,  AT 
SA/H  =  ioAlAT 
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c:  =  siK  (0.s*alo^) 

ij  =  0"r 

Afl7  =  WUATH-PI ATA 
CAXH  =  (  1  .G-0)  *SU'  (  AA7) 
i)  =  (SMtt) 

CA/A  r  CA,/M«0 
CA/h  -  CA/h-O 

r  w  r  ( ‘-.A7AeCA7A  ♦CA/A-CA/A) 

::  f^LAToASIf  (.^) 

SA/A  =  SA7A/0 
CA?A  =  CA/A/0 
SA7b  =  SA/H/O 
CA/H  =  CA7H/0 
HA7A  =  ATAK'2(SA7A.CA7A) 

,VA/e  =  ata^2  i<;a/h  .CA/H) 

ALON  =  IVA7H-HA7A 
If  (A(-5  (  ALON) -P  I )  120»10S«105 

lOS  IF  (AION)  110,110.115 
110  PA7H  =  PA7R*TW0F“I 
(.0  TO  1?0 

1  15  PA7A  =  H47A»TI*PPI 
120  hlAT  =  CF‘^Q*E50*AI  AT*ALAT 
HLON  =  AO/SORT (HLAT) 

HLAT  =  HLnN«CF50/HLAT 

HLON  =  HLO^'^ALAT 

CLAT  =  CA7A/HI AT 

CLON  =  SA7A/HLON 

HLAT  =  CESQ*E5C*HLAT*BLAT 

HLON  =  An/SORT(HLAT) 

HLAT  =  HL0N*CFSU/HLAT 
HLON  =  HL0N*8LAT 
ALAT  =  CA7H/HLAT 

HLAT  =  ( 3.n*AA7/5Me-ALAT-?.0»CLAT) /SHH 
Al Of  =  SA7H/HLON 

BLON  =  ( J.O»BLON/5MB-ALON-?.0*CLON)/SHb 
AA7  =  3.0*SHH 

ALAT  =  ( ( ALAT-CLAT) /SHB-2.0*HLAT)/AA7 

ALON  =  ( (ALON-CLOM) /SMe-2.0*BLON)/AA7 
AA7  =  (RA7R-RA7A) /5MH 
ENTRY  PCnORD 

rLATP  =  (  (ALAT*5MP*t<LAT)*SMP*CLAT)  *5mP*RLATA 

HLONP  =  ( (ALON*SHP*HLON)*SHP*CLON)»SMF»ALONR 

WA/P  =  AA7*SMP»Ra7A 

PFTURN 

END 
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•DECK  GETELV 

SUBROUTINE  GETELV (L AT  IN .LONIN . AMP ♦f A? .FLVTN ) 

C  PURPOSE  * 

C  COMPUTES  ALTITUDE  AND  IMPEDANCE  AT  A  GIVEN  LATITUDE  AND  LONGITUDE. 

C  ASSUMES  A  DATABASE  TAPE  WITH  ?40  RECORDS  PER  /DNE  OF  POINTS.  60 
C  WORDS  PER  RECORD.  AND  2  DATA  POINTS  PER  WORD.  IhesE  GROUPINuS  COVER 
C  THE  FOLLOWING  MAP  SECTIONS* 

C  A  zone  SPANS  1  DEGREE  LATITUDE  (SOUTH  TO  NORTH)  AND  ?  DEGREFrS 
C  LONGITUTDE  (WEST  TO  EAST)  OF  THE  MAP. 

C  A  RECORD  SPANS  1  DEG.  LAT.  (S.  TO  N.)  AND  30  StC.  LON.  (W.  TO  F.) 

C  OF  A  ZONE. 

C  A  WORD  spans  I  MIN.  LAT.  <S.  TO  N.)  AND  30  SpC.  lON.  (W.  TO  F.). 

C  OF  A  RECORD. 

C  A  POINT  LIES  ON  EVERY  30  SEC.  LAT.  AND  EVERY  30  SEC.  LON. 

C  ALTITUDE  IS  INTERPOLATED  WITHIn  4  POINTS  WHICH  LIE  NORTHEAST. 

C  SOUTHEAST.  NORTHWEST.  SOUTHWEST  OF  THE  INPUT  POINT. 

C  inputs: 

C  LATIN  =  latitude  IN  DDMMSS  FORMAT. 

C  LONIN  =  LONGITUDE  IN  ODDMMSS  FORMAT. 

C  outputs: 

C  AMP  =  IMPEDANCE  AMPLITUDE. 

C  FAZ  =  IMPEDANCE  PHASE  (RADIAN). 

C  ELVTN  =  altitude  (METER). 

c  calls  routine  index  fori 

C  INDEX  s  POINTER  TO  A  RECORD  IN  THE  DATABASE. 

C  IXPT  =  POINTER  TO  A  HALF-WORD  IN  A  RECORD  IN  THE  DATABASE. 

COMMON  /MAP/LATSW.LATNE.LONSW.LONNE 
COMMON  /GETELV/L(120)  .IMPdZO) 

COMMON  /UNPACK/NPACK (60) 

COMMON  /FLAGS/IXPT.NFLAG.NRPC.NRPDB 
COMMON  LAT(99R) .L0N(999) .MSINOX (5800) 
data  KPASS1.MASK8.MASK11  /I .3776. 3777B/ 
data  INDDIM  /5800/ 

C  READ  and  store  DATABASE  TAPE  IN  FIRST  PASS  ONLY. 

IF  (KPASSl  .NE.  1)  GO  TO  630 
KPASSlsO 

NRPC=240* ( (LATNE-LATSW) /lOOOO) 

NRPD8=NRPC* ( (LONNE-LONSW) /20000 ) 

IF  (NRPDB*!  .LE.  INDDIM/  GO  TO  610 
PRINT  9700.  NRPDB  5STOP 
610  CALL  0PENMS(2.MSINDX. INDDIM, 0) 

DO  629  IRT=1. NRPDB 

BUFFERINI 13.1) (NPACK (1) .NPACK (60) ) f IF (UNl T ( 1 3 ) >620.629.620 
620  CALL  WRITMS(2. NPACK. 60. IRT) 

629  CONTINUE 

c  IF  latitude  and  longitude  lie  outside  map.  previous  data  are  used. 

630  IF(LATIN.LT.LAT5W  .OR.  LATIN. GE .LATNE )  GO  lO  190 
IF (LONIN. LT.LONSW  .OR.  LONIN. GE .LONNE )  GO  TO  190 

C  CONVERT  TO  RELATIVE  LATITUDE  AND  LONGITUDE. 

nflag=o 

LATX=LATIN-LATSW 
IF  (LATIN. GE.O)  GO  TO  200 
C  CORRECT  FOR  NEGATIVE  LATITUDES. 

IF  (M0D(LATX.  10000) .EG. 0)  Go  TO  200 
LATX=LATX-4000 

IF  (M0D(LATX. 100) .EO.O)  GO  TO  200 
LATXsLATX-40 


34 


t  Mf/SITIVf  (.ATlTUOfcS. 

fUa  Lfll  St  (-  =  ■^00  (LATX  .  1  On  ) 

L  AT(>w  =  t  A  TX  -L  ATSt  (. 

IF  (L  TFFf  .i.;F  .  30)  C-O  TO  t'AO 
(_ATS.:i  iT'1^■ 
a  Tf.'-i.  _  F'ly  »  30 
b '  ■  T '  ’  .  ‘ 0 

'  .)  .1  T  =  i  ;  ;M.  3fi 

a  7  iOo 

[  A  (  L  6  TDA-  .  1  Ff-r.  li ) 

IF  (LAT^.F^j.S'vOO)  L  ATN=LATDM*i.lOO 
ehO  LOFiX  =  Lr)MN-LON'S« 

IF  (LOMlN.Gt-.O)  bO  TO  300 
C  CCiB^^tCT  FOP  NEGATIVt  LONGITUDES. 

IF  (MOP  (LOt.'X  ,  1  0000)  .FO.O)  GO  TO  300 
LONX  =  Lr'NX-4000 

IF  (MOD (LONX  ,  1  00)  .F 0.0)  GO  TO  300 
LONX=LONX-40 

C  OFfOCEFt)  WI7M  POSITIVE  LONGITUDES. 

300  LOr'SEC  =  MOP(LOF.X.  100) 

LONDm  =  l  ONX-LOF  SEC 
IE  (LONSEC.GE.30)  GO  TO  340 
LON*(  =  LONOM 

LONE=LONDM*30 
GO  TO  3b0 

340  LON(xi=LONDM*30 
LONE=LONOM*100 
LONM=MOO(LONDM, 10000) 

IF  (LONM.EO.SSOO)  L0NE=UONnM*4100 
360  continue 

C  GET  A  DATA  POINT  FiO*»THEAST  OF  THE  INPUT  POINT. 

INDNE=INOFX (LATN.LONE) 

IF (NFLAG.EO.n  GO  TO  120 
CALL  WEADMS (?.NPACK.60« INONE ) 
call  unpack 
INF=L (IXPT) 

C  GET  3  POINTS  SOUTHEAST.  NORTHWEST.  AND  SOUTHwEST  OF  THfc  INPUT  POINT. 

IE (MOO (LATN. 10000) . EO.O. OR. MOD ILONF ,20000) . tu . 0 ) GOTOl 20 
C  ALL  ROINTS  lie  IN  THE  SAME  /ONE. 

ISE=L(IXPT-1) 

CALL  PFAONS (2.NPACK ,60, lNDNF-1 ) 
call  UNPACK 
INW=L (IXPT) 

ISW=L ( I XPT-1  ) 

GO  TO  1«0 

C  THE  4  POINTS  lie  in  DIFEERENT  /ONES. 

120  CONTINUE 

INDNW=IN0FX (LATN.LONW) 

IF (NFLAG.EO. 1 )  GO  TO  140 
call  PFADM5(2.MPACk,60,INONW) 
call  unpack 
INW=L (IXPT ) 

140  CONTINUE 

IN(;SE  =  TNDf  X  (L ATS, LONE  ) 

IF  (NFLAG.FC  .  1  )  I’O  TO  IKO 
call  re  AOMS  (2.^.PACK  ,60  ,  INI.SE  ) 

call  unpack 


ISE=L ( I*PT) 

160  CONTINUE 

lNf)SW=INDEX  (LATS.LONW) 

IF (NFLAG.EO.l)  GO  TO  100 
CALL  RE  Ar)MS(2,NPACK,60,lNO<;w) 

CALL  UNPACK 
ISW=L ( IXPT) 

160  CONTINUE 

AMP=FLOAT (SHIFT ( IMP ( I XPT ) «-l 1 ) . A. MASKS ) * . 00 1 
FAZ=FL0AT (IMP (IXPT) .A.MASKI 1) *.001 
IF (NFLAG.FO.l )  PRINT  9702 

ELF  AST  =  FLOAT ( ISE ) ♦FLOAT ( ( INE-ISE) *MOO (L ATSEC . 30 ) ) /30  . 
ELw£ST=FLOAT(ISR) ♦FLOAT <( INw-ISW) *MOn (L ATSEC t 30 )) /30. 

ELVTNatLWFST^ ( F LE AST-ELWEST ) *FLOAT ( MOD ( LONSt C , 30 ) ) /30 . 

RETURN 

190  PRINT  9701.  LATIN, LONiN 
RETURF' 

9700  FORMAT  (*  DATABASE  HAS  *110*  RECORDS.  INCREASE  ARRAY  MSINOX  AND  VA 
IPIABLF  INDOIM  TO  ACCOMMODATE.  RUN  ABORTS.*) 

9701  FOPMAT(*  *,4H****,*LAT=*I10*,LON«*I10*  REOuESTED  ARE  OUTSIDE  MAP  R 
lEGION  -  NO  elevation  RETURNED*.//) 

9702  FORMAT)*  . . INDEX  REQUESTED  EXCEEDS  SCANLINES  GENERATED.  TH 

lUS  PREVIOUS  ELEVATIONS  WILL  PE  USED*) 

END 


Index 

Function  index (lat.lon) 
c  PUPhOSE 5 

L  computes  pointers  into  DATAPASE  for  a  given  laTITUOF  and  longitudf, 
C  IF, puts: 

C  LAT  =  REIATIVF  LATITi.CF  In  ODmmsS  KODuat. 

C  LON  =  RFlATIVF  LONGITUOF  IN  ODOumSS  FORMAT, 

C  NRPC  =  NUMBER  OF  RFCORUS  PFR  COLI'mu  OF  mAR  APFA  2  f'F.GPEES  LONGITUDE 
C  WUiE. 

C  x,PPi)P  =  NUMBER  OF  wF.  roPOS  IN  THF  (JATAhASE. 
c  Outputs: 

C  lF.Utx  =  POINTFP  TO  PFCOPO  IN  DATa-ASt. 

C  IXPT  =  POINTER  TO  (lUft  POINT  (=  HALF-wURU)  In  QATAPASE  RECORD. 

C  f.elAG  =  1  INUICATFS  computed  RECOPO  POINTFP  IS  OUT  OF  RANGE. 

C0MMl)^:  /slags/1  xpt.nflag.nprc.nppoh 

LATS=no'')(LAT.1C0) 

LATm  =  m(V)  (L  at . 1 ooon ) -LATS 
LATf)=  U.AT-LATm-L«TS)  /  looon 
LATM=l ATM/lOO 
LUUS=M()D  (LON.  100) 

L(JF  M  =  M00  (LON.  1  0  000  )  -LOF'S 
L('NO=  (LON-LONM-L  Df.S)  /  lOOlTO 
l(ir  rivjM/l  00 

11  OM)=  (LOr  U/2)  *UPPr  .MOf'  (L('f  D.2)  *l2').i  ONm* P. LDNS /■  3(' ♦  1 

IM'E  X^lLO:  ^-'♦LATu*?P0 

I*FT=L«TM«2.LATS/3n.l 

IF  ( IFP  F  X.GT.NPPUP)FFLAG=1 

PFTUPF 

F  F.l) 
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"(iECr  AC* 

;T  TNJI-  i  ... 

C(1I.',V0''.  /  1  U  \/  /  L  !  i  r'-'’  )  A  !  ‘-'p  (  1  r' '•  ■■ 

/I  :\;.aCK /N- '  •<  (.‘.'i; 

,  .'TA  ‘.-a  1  1/ j7f  1-^/ ,  l'J/1  /  y  ;j  tlrK/ 

I,  (  ,  a:'  I  .-Ift- 1  ),  I  i,  ,  V4,  , 

L  ■:  >■  ;  -  w'':i '  s>‘  i  ;  <  Ah,  1  I-  T  ( Me  ,  -l.  A  I  i  •>  r 

j  ‘.ip  i  M  )  =  AM'  1  AAsk  1  ,  .  ^  -I  I  P  T  (  'jP^  ‘.'A  ,  _<.••  I  ) 
t_  ( *.M  !  ;  =  af j( :  ( M A  5jK  1  ;  •  Ah, T  P  f  { !.,!-'»  -■  A  ,  -  i  >4 ;  )  o  >- 
I  K'p  (  M  »  ',  )  =  AM.'  (  w  A  AK  1  S  ,  •'il'.K  ^  A  1 
f'Ci  CriNTlA'jP 
►yfc  T'jSf. 
tM' 


^OECK  IKT 

SUBROUTINE  INT  (I*(‘| 

C  inputs: 

C  I  =  position  IN  X  AnP  V  awPAYS  on  which  to  center  CALCULATIONC 
c  K  =  position  in  aprays  ;  and  /p  to  stork  calculatfp  values 

COPMON  /()RNDl/7PP.XO(E)  .7  ( lOOv)  ,7PD(<S)  ./P(1009)tXD(S).KIlr'nR) 
IM0=I-1  ■tXIMOsX  (  IPO)  S7IPC=7(1M0)  A7!  =  /(1) 

IP0  =  I»1  SXIPO=X(IPO)  $Xl  =  X(I)  «txK  =  X(K) 

C=( (? (IPn)-/lM0)-(/I-7IM0)* (X]DO-XlMO)/(XI-XlPO) ) / 

1  (  (XlPC»XlPO-XlPO»XIMO)“(Xl»XI-XlMOoxlPO)* 

2  (XlPO-XlMO)/(XI-XlMO) ) 

b= (  (7I-7IPO) -C* (XI«XI-XIV0*XIP0) ) /<x I-<]MO) 

A=7I-XI*(B*C*XI ) 

Z(K)=A*XK*(6*C*XK) 

7P(K)=R»2.*C*X)< 

ZPP=2.*C 

RETURN 

END 
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<s"r--,in_F  -> 

If  rrpF 

SUPCfOdT  T^^  I'.tO/'F 

r  oFFFoFNrF<;; 

r  (1)  J.t- .  .  C.<>kp>- C  T  MMS  fryi- ■/ 

r  If 'hO^OC’e  fir  I  iUS.  I  IWf  TFi-PAIr..,  FsC,'-  >-tr'''FT 

r  TFP  'iT-nsi.  ‘->h,  'juvf-MKP.  lSlfc7. 

COMMOr.'  /PlS/T/,dPI  ,PI  .t-iAFPI  ,'IPTPI 

CO“t.‘OM  /OPf  I'  1/  /pp  .  (h  >  .  /  (  ,  /  t  11  (H  .,'P  (  1  Ou'^  )  ,  (  b  )  t  <  (  1  0  Il'-I ) 

COPI^CA  /(..PM-V/pO  .Orrj.K.ij.O’MX  ,P<?,|'f  (.1 
CCMMor:  /Tr //ppi  ( lOus) 

COMMON  /70TA/XPPAY ( lb) 

common  /SfiPUI  /  S  (bPp)  .OP  (PPP)  .01  (MSP)  ,LLM 

COMMON  /CnCT  /  1  TCT»LT0P.r'PTC1  .  ?Xmt P  .  7m  I N .  AL T  bv.i  (  3 )  .  TSw  (  .) 
r  70.  7Pn  ANO  XD  PP'JVU't  STCPAOI  fop  rjfOAlIVF  XPi^. 

oIMFNFTON  T/tP  (  1  )  .  T  (  1  000  )  .<i  »  (b)  tOX  (S)  .r!  .1  (9)  .r  X  (  p)  .  mu  (b)  .XS  (  p) 

CiIMrNfTo'i  »!F  (9)  .<M  1000) 

COMPLF  X  OF  L  I  »OtN.r  LK  AF  .1  1  .r  ?.F  3.f-l  ,  <  V  . '’3  .  SIIm  .  T  t  . 

1  w .  wr .  P 0  .  w  1.  4  X  .  T  b  .  r  I  N( I  . f  I 

PATA  KPAbSl/l/ 

OATA  ((0X(N),  n  =  1,  b)  =  .«)^'7163br^br^l  ,  .Ib/wsissyb.  .i.blb'' ?3f- 1  b, 

1 .7A83346?t>4.  .948b93976F’) 

OATA  ((Ow(t<).  K  =  1,  b)  =  .b9104P4494.  .S3rib3i4>bb  * 

1  .?9ePn2NS«rt3.  .  1333426bPb) 

OATA  ((FX(t),  n  =  1,  ..)  =  .10bl4(i2a26.  .i/fc^/Ablbb*  ,  bspi-b  -  o )  2b . 

1  .9373342493) 

OATA  ((FPlf).  ■(  =  1,  b)  =  .0f''3f>«0bl9b9.  .19Ku9b2bbb.  .  2'^  ?b  2  7  34Sb . 
1.I9?36’03S7) 

CATA  (T7f.P(l  )  =  l,0) 

9400  FORMAT  (<*  AlKbOPNt  TARGET  ALTITUDE  (4r9.2*)  l.ltS  -iElOv  TEPb£.lN  klE 
IVAirnr'  (<rE9.c").  Mi/N  AHOPTS.®) 

FTlWAVFb)  =  CmPLX(C(IS(9AvEP)  «  SI  M  (  a- VE  P  )  ) 
r-----------------  ------------------ 

r  START  first-pass  INITIALI/ATIDN. 

IF  (KPAcSI  .f.E  .  1)  GO  TO  9999 
KPASSl=n 

C  FPf:;;  IS  licc-r,  If,  ,Vf.^  JK.  This  PPOOPAM, 

FRFC  =  ARRAY  (  1  )  <*1  ,E -3  YRA|j  =  APRAY(13)  SFlAT  =  ARpAY  (  1  1  ) 

FORM  =  array  (iH)  y-SmOOTH  =  ARRAYdb)  IFTA  =  ARpAYI?) 

fiTFPOs  9 

T  (ri7F09)  =1  .0 

r  ,i.AvF  r.tlx'PFR  IN  ATMOSRRFPt  AT  SOREaCE  OF  EARTh 
<lAVF  =  ARRAY  (  1  (I  ) 

AMICPO  =  1 .0/ ( TwOPI*FpEQ) 

TX  =  S9PT  (E -t  i-bFTA  )  ".040R3b9S49 
r  If  lTIALI-’F  ''I'ANTlTltS 
X  (  1  )  =  0  . 

ii.  (  1  )  =  1  . 

C  Tt-AY  ciYOOTp  p(,|,N.,.  FAbli-  IF  SmO0TH=1. 

IF  (SMOOTH  .ri,.  1.)  7(  1  )=7P(  1  )  =ZRE  1  (  1  )  =0. 

oor.^  CO'ITIN'.lF 

o  FnO  EIPST-paSS  IIn]  riALi/ATlON. 


r  STAE  T  INITI  ".L  I/ATION  E  uP  CURPr  NT  PATH. 

r  o.-mri'TF  DIOTAfCr  A-li.i  OloIaNCE  INI'pEMEmI  (hm)  TD  TAxolT. 
ri-Avix  (f.PTs)  /looo. 
r  ixic  =  x  (?)/iouri. 
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LTOP  =  0 
MOI-lAT  =  1 

C  iMPfcOANCf  'SHOULD  PtFFXHNCKti  To  100  Kh/ 

L  KPFOUlNCY  UFOFNOENCF  01  IMPEDANCE  =  *SDPT ( F PFC / 0 . 1 ) 
otLi  =  CMPLX(np(i),  01  (in*«.oPT(FPEo/o.i) 

IMOST  =  DMAX/OINC  ♦  .01 
ulLTX  =  1000.*nlNC 
T(l)  =  1 ./SOFT (Dfc L TX) 

C  FLAT  MPTH  Theory  OtSCPHED  IN  PFFEREMCt  3 

on  Ji"  K  =  1  «  E 

L  yTTr  rguA  UON  f  OP  FLAT  EAPTH  THEORY 
X(-K)  =  OFLTX"OX(K) 

call  OPOUND  (  -K.  ?.  0.  0) 

3?  l«iG(K)  =  FLEAF(^AVEg  0..  0,.  X(-K),  OF  L  I  ) 

1  =  2 

C  END  IMTIALI7ATION  FOP  CURRENT  PATH. 

- 

C  START  LOOP  ON  SURFACE  POINTS  ALONG  miiFFOhD  HASELINE. 

37  LOMTlt'UE 

C  TREAT  SMOOTH  POUND  EARTH  IF  SmOOTH=1. 

IF  (SMOOTH  .EO.  1.)  /{n=/P(n=7PPl  (I)  =  Ti. 

7PP=7P01 ( I ) 

T(I)  =  l./SnRT(X(l)  ♦  OELTX) 

OMX  =  nxicos  (  X  ( 1 ) /PAD) 

I  GO  -  ?  -  MOD (I,  2) 

IL  =  I  -  100  -  1 
C  EOUATION  2.1H,  REFERENCE  1 

RO  =  SORT (?.*RAU* (PAD  ♦  2(I))*OMX  ♦  7(1)**?) 

IF  (NOFLAT  .EO.  2)  GO  TO  AS 

c  Flat  oetefmines  distance  to  which  flat  farth  theory  is  used  rffore 
C  irregular,  InhOMOGFNEOUS  terrain  )AAY  be  inserted. 

AO  IF (  1  .lF .A.OP.FLAT.GT.X ( I ) )  AI.aS 
A1  w(I)  =  FLEAF(WAvf,  0.,  0.,  «(I).  DELI) 

GO  TO  PO 

L  SPHERICAL  FAPTH  THEORY 
AS  SUt-’  =  0. 
noflat  =  2 

DO  so  K  =  1 ,  S 

C  R,  DEN  and  II  ARE  FROM  f, POUND 

call  C-F0"(I,  -K.  1,  0.) 

C  EcAJATION  ?.27,  REFFRFf'CF  1 

term  =  wG  (K  )  *F  I  ( -.’A  VE*R  )  *r)FN 
so  SUM  =  (U*Gw (K) ) /S' RT ( X ( I )  -  X(-K))*TFRM  ♦  SUM 
C  SIMPSONS  PULE 

SUM  =  3."T(1)*SUM 
KK  =  1 

IF  (IL  .LT.  3)  GO  TO  100 

DO  60  K  =  3,  IL 

CALL  GFOM(l,  K.  1,  0.) 

TF.PM  =  u*T  (IX-1  )  *»i  (N  )  *EI  <-.'AvF*P)*DEN 
IF  (Kw  .EO.  2)  GO  TO  SS 
S(JM  =  A  .*T  (  I-K  )  *TF  PM  ♦  SUM 
IXK  =  ? 

GO  TO  60 

SS  SUM  =  ?.*T ( I-K ) *TF PM  .  SUM 
KK  =  1 

h('  CONTlruP 
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i  100  continue 

'  C*LL  OfcOMd,  ?,  1.  n.) 

SUV  =  T  ( l-<;  I  *T  (  1 1  (?)  *E1  (-1'.AVE*V)  *l't *  S'.IV 

call  GE0M(I.  I  -  If-0,  1.  0.) 

F?  =  U«T  (  IL  )  *W  (I-U’O) ‘t  I(-V<A\/£*R)  "UEN 
SUM  =  (Sum  ♦  T  (  1 00)  *F? )  * .  3  31333:)»DF  LT  ». 

IF  (ICn  .£(,.  ?)  (>n  TO  hh 
El  =  E? 

F?  =  TFRM 
GO  TO  70 

6ft  call  GtOM  ( 1  t  1  -  1 .  1  ,  0.  ) 

El  =  li«T  (  r-2)  ( I-l  »  I  (-sA^F 

Sum  =  SUM  ♦  .0AJ33333333»OFlTx‘ (S.*T  (  1  >  "E  1  ♦  t.*T(?)«-E? 

1  -  T(3)*TFRM) 

70  0  =  TX/T{ 1-1) 

RHO  =  1 .  ♦  2P ( I ) *o? 

RHO  =  7PP/EHO 

TERM  =  1  .?/T  ( 1 )  «T  (I-l  )  »LMPl  X  ( 'j.O  “  OF  LI  ♦  Cmdlx(0.«  -wO/wAVE)) 
WX  =  1,  -  CMPLX (0.0)* (SUV  ♦  ?,/T ( 1 ) * ( .Afeftftftftftft7*F  1 
1  -  .066ft6ftftft7*E?) )/ ( 1 ,  .  term) 

l«  ( 1  )  =  «XX 

90  OIST  =  .001*X(1) 

CALL  INOF  (0. .X ( 1 ) .EOPH.FIND) 

TS=W ( I ) *EIN0 
PhA=CANG( TS) 

C  PHIC  IS  Tvt  SECONDApy  pmaae  CT'A^fctIon  in  paijIAnS. 

PHIC  S  -  (PHA  -  yAVE»(R0  -  xd))) 

C  sec  is  the  secondary  phase  CORRECTION  Ik  V I CpOSFCONPS . 

SEC  =  PHIC*AMICP0 

C  SAVE  SFCONDAPY  phase  factor  at  TAwoE  T  (’POUND  P0l\'T. 

IE  d  .NE  .NPTCT  )  Gt)  TO  9h7 
TSiAdTCT)=S£C 
R97  CONTINUE 

I  =  I  ‘  1 

IE  d. I  E. OPTS)  .,0  TO  i7 

C  EnO  loop  of;  SUPEATE  POPjTR  Ali.No  h:iFF')Pi)  HASr.LlNfc. 

C-------------------- . 

C  START  AlPHORNE  TAP;-' T  r  T-p, ,  t  A  T  !  . 

IE  (2F'IN.LF  ,0,  ) 

LTDP  =  1 
E?=(0..0.) 

G?~ ( V . . 0 . ) 

OLDR=0. 

II  =  I  =  NPTCT 

IGO  =  ?  -  VOC'd’  ?) 

00  iSS  <  -  if  ^ 

XdOOO  *  K)  =  xd)  -  OF  L  T  j'T  >.  (-v  ) 

XS()')  =  «(100  0  *  •') 

1S5  call  GPOUNO  '  lOOD.i*  .  II  .  f'. '■  ) 

call  ( rEV''F:N  (  «  ,  M,  l'-n)ST,  XS.  -E.  “  ’  S.  1) 

HL  =  /NIN  -  /XMTP 

L  APORI  IE  TARGET  ALT'dUOF  TS  EF  I  'i"  TERRAIN. 

IE  (HL  ,ST.  7  (I  )  ;  GO  TO  I'-A 
PRINT  t,.,on  .HL  ./ d  )  SSTOP 
1S4  7/  -  PAO  ♦  HL 

RO  =  SORT  (  ^  .  »R  AT’*:- 7»TN  C'Ns  (  *  (  I  j /RAO  )  .  hl**,-) 

SliM  =  T:  , 

i 


\ 


■)0 


lJ()  1S7  =  1  • 

CALL  i)vi  {  1  ,  _n  .  V  .  -tL  ) 

Tew^■  =  A,!-,  (r  )  <-f  I  (  -  •  I ' 

1S7  SIN  =  laS'.-  1“  )  -Tr -'‘•■/S''»aT  (X  (  1  )  -  >(-:■.))  ♦  Si'' 

C  SIMPSf•^S  i^liur 


SUi"’ 

=  3.<» 

T  (  1  )  *Si .  •> 

C&IL  bKnf/ 

{  1  ,  1  ,  8  ,  HL  ) 

FI 

=  T  (  1  ) 

«T  (  I  -?  >  Cll*')'^  ’J*1A  (  ?) 

G1 

=  F  I  (- 

4  V  L  *  P  ) 

SUM 

r  c  1 J  M 

*  F  1  *(a] 

U  0 

=  1  - 

luu 

K(-0 

=  1 

IF 

(  I  L  0  . 

LT,  3)  GO  TO  1  IS 

on 

18?  K 

=  3,  ILO 

F  ;< 

=  F? 

o3 

=  G? 

F? 

=  FI 

<■.? 

=  G1 

CALL  GFO«(I.  K,  h),  ) 

G1 

=  FI  (- 

■vA  Vt*P) 

FI 

=  T(K- 

1  )  *T  (  I-K  )  *U*[IFN*.>I  C*  ) 

(lELTO  =  .■JAVt*(P  -  OLOP) 

IF 

(KGO-?)  158,  MS.  17? 

158 

Sum 

=  SlJM 

♦  4,*F1o(a1 

KUO 

=  ? 

GO 

TO  18? 

150 

SUM 

=  SUM 

*  ?,*fi*gi 

KGO 

=  1 

IF 

(APS(OElTG)  .GT.  .?)  170,  18? 

170 

KOO 

=  3 

SUM 

s  SUM 

-  F1*G1 

GO 

TO  18? 

C  EOUATKiN!  ?.3S,  wtfFPEN'CF  1 

17?  SUM  =  S!!"*  *  CMPLX(0..  3./0tLTu)  *  (Gl*  (F  1  -  .  S*  ( F  1  *C  N'PL  *  (  2 
U.  J.*OElTG)  -  A.oF  ?»CN'PLX  ( 1 ..  OELTG)  ♦  F3*CmPL>(2..  ()fLTr< 
?)  )/r'ELTG»*?)  -  6?«(F?  -  .S*(F1*CMPLM?.  .  UFLTC  )  -  A.o 
3F?  ♦  F3«CMPLX(?.,  -OELTG) ) /UELTG"*?) ) 

18?  OLOH  =  P 
US  CONTINUE 

IF  (IiaC  .FO.  1)  GO  TO  iFf 
F3  =  F? 

G3  =  C-? 

F?  =  F  1 
G?  =  G1 

CALL  GFOM(I,  1  -  1.  ?,  HL) 

Gl  =  F  I  (-'a4VE*P) 

FI  =  T  (  I -?)  *T  ( 1  )  *U<^0F,N**»  (  1-1  ) 

OFLTG  =  WAVE*(P  -  OLOP) 

IF  (AfS(OeiTG)  .UT.  .?)  1H4.  MS 
C  Ef.L'ATION  ?.3S.  PKFPPtt'LE  1 

1(J4  SUM  =  SUM  ♦  CMPLx((i,«  l./Of  LTG)<»(G1*<F  1  -  .S*  (F  MCmplx  (? 

1..  3.»OFLT(-)  -  -‘.n-F  ?4CMPL>  <  1  .  .  DELTO)  ♦  F  ,T»Cmpl  >  t  c  •  •  OFLI'- 
?)  )/UELTG*"7)  -  iit'^CF?  -  (F  MCMPLX  (?.  .  fFLTC)  -  A.« 

3F?  ♦  F?*CMPLX(?..  -ncLTG) ) /OFLTG**?) ) 

GO  TO  ISM 

jtju  5^.^<  =  tljM  ♦  l.?S*F!*Gl  ♦  ?.*F?*G?  -  ,?S*N  J*G.3 
188  DO  100  K  =  1.  4 
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C6LL  Grop'd,  loon  ♦  k,  3,  hd 

9|iM  =  SUM  ♦  1,  (M  *Vtr' (M  *P  1  (-Wftvf  “Pi  »r)F.N*U/ (  S‘J'»T  (  XS  (K  )  )  *  1  (  1  )  ) 

n  GG-  jT  I  Ml  ir 

=  ,3333.1  )  lT  x«SUM 

G  =  TX/T(I  -  i) 

i<.v  =  (1.  -  C-v-LX  Id.'.)  «SUM)  *.S 

C^Ll  1'  (Ml.x  (  1  )  .tOkH.X  IM'D 

F='A''  =  -  (l.AMu(TS)  -  MAVt*(W0  -  X(l))) 
r  r-<<ro  IF  jiJr  IM  MICP06ttO‘-/OS  wHltr  F»AD  IS  THc  PHASf  IN  RADIANS 

fr.-SGC  =  F  PAl"»A*xICPO 

r  ti/F  phAcf  at  lApioFT  POINT  F  Ow  TIME  niFFFPENCF  CALCULATION, 

.'LTSV  < ’TCI  )  =  F  XStC 

PFTI'SM 

(  f  p 


*OtCK  GROUND 

subroutine  GROIINUd*  K»  IGO,  MH) 

C  THIS  SUBROUTINE  USES  VARIABLE  IMPEDANCE:  AND  TERRAIN  TO  CALCULATE 
C  GROUND  WAVE  FIELDS 

COMMON  /&RND1/7BP./U  (M  .7(1009)  »/RO(f>)  . /P  ( 1  009 )  ,  XD  (B  )  .  X  (1009) 
Common  /GRND2/(-0.t  I-  N.p.U.0MX.R2«0ELTA 
COMMON  /70TA/APRAY ( IS) 

COMMON  /FOROI/S (999) .DR (999) .01 (999) .LLM 
complex  delta, pen 
c 

call  INT  (K.I) 

FRED  =  array ( 1 ) *1 .E -3  TnAVF  =  ARRAY(IO)  TA  =  ARRAY(13) 

RETURN 

C 

ENTRY  GEOM 


HIT 

=  7(1)  ♦ 

HH 

IF 

( 

I  .FQ.  K) 

GO  TO 

T  = 

(  X  (  I  )  -  X 

(K ) ) /A 

GS 

s 

A  ,  7('<) 

GX 

= 

A  »  HIT 

CT 

= 

COS(T) 

ST 

= 

SIN(T) 

OT 

= 

OMCOS(T) 

C  EOUATIDN  ?.?3,  REFFREE'CE  1 

R2  =  SORT (?.*6S*GX*nT  ♦  (HIT  -  7(K))*»2) 

c  eouation  ?.22.  reference  1 

C  E(JUATinN  2.30,  reference  1 

R1  =  S0RT(2.*A*GS*0MC0S(X(K)/A)  ♦  /(K)<»*<f) 

R  =  R1  ♦  p?  -  po 
C  EOUATION  2.30.  REFERENCE  1 

C  U  =  PROJECTION  FACTOR  TO  THE  SURFACE  OF  If TEGrA T I  ON . 

U  =  X(«)*RO*SOPT(1.  ♦  7P(K)»*2)/(Rl*P2*x  1  I)  ) 

IF  (IGO  .LT.  3)  U  =  U*(*(l)  -  X(K)) 

C  EOUATION  2.30,  REFERENCE  1 

C  PD  =  normal  partial  derivative  to  the  surface  OF  The  Ijii-OUNO. 
PD  =  (A*nT  ♦  7(1')  -  H1T*CT  ♦  Gx*7P(K)  6A*ST/f>S)/R2 
XK  s  X (E ) 


IF 

(XK  .GE. 

S (LLM) ) 

GO 

TO  12 

LLL 

=  1 

IF 

(XK  .LT. 

S  ( 1  ) )  GO 

TO 

10 

LLMO  =  LLM 

- 

1 

IF 

(LLMO  .L 

T 

.  1 )  c-a 

TO 

100 

DO  13  LLL  =  1.  LLMO 

IF  (XK  .LT.  S(ILL  .  1)  .AND.  XK  .ON.  S(LLL))  GO  TO  10 
13  continue 
100  CONTINUE 
GO  TO  21 

C  frequency  DEPENDENCE  OF  ImpEiiA»^CE  =  *SoRT  (FRE<)/0 . 1) 

C  DELTA  =  COMPLEX  LOCAL  (^ROUNU  IMPEDANCE  THAT  IS  ADIjED  TO  A 
C  FUNCTION  INVOLVING  wAVE .  R2  AND  PO. 

10  UELTA  =  CMPLX  (OR(l  LL)  ,  U I  (LLL  )  )  *Sl)RT  (F  Rti:  70 . 1) 

GO  TO  21 

C  FOP  VALUES  OF  x  GPEATEP  THAN  THE  LAST  S  VALUE  READ  IN,  THt  UA<;t 
C  IMPEDANCE  IS  USED 

C  FREQUENCY  DEPENDENCE  OF  IMPEDANCE  =  *SORT (F RE M/n . I ) 

1?  UELTA  =  CMPLX  (DR(ILM)  ,  DI  (U  «)  )  *Sc.wT  (FRE'-'/O  .  1  ) 

21  CONTINUE 


C  delta?  =  OEPItfflTIVF  PAkT  OF  U'PEOANCE 
C  tfiUATICN  ?l,  PEftt'F''lCfc  2 

delta?  =  CMHLXd..  -1  ./ (wAVE*E?)  )  *P(, 
C  tOUATIOi'*  ?0.  <iFFEKFNCt  ? 

IS  DEN  =  OFLTA?  *  DELTA 
weTUf<Ei 
^0  IJ  =  0. 

W?  =  (,x  -  A  -  7  ( I  ) 

Pt  TUWE' 

END 


*UECK  OMCOS 

FlINCTIOO  OE!COS(A) 

C  OMOCOS(X)=  I  -  CD'S)*) 

IF  (Ahs(X)  »GT.  .lb)  GO  TO  40 
IF  <X  .EO.  0.)  GO  TO  bO 
S  =  X*x 

0“C0b  =  T  =  .b*S 

B  =  4, 

10  T  =  -T*S/r»*(B  -  ].>) 

OMCOS  =  OMCOS  ♦  T 

IF  (T/OMCCS  .LE  .  .^E-S)  GO  ti'  SI 

w  =  B  ♦  2. 

GO  TO  in 

40  OMCOS  =  I.  -  cns(X) 

BETUBEi 

50  OMCOS  =  0. 

51  RFTUBM 
EM) 
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“UtCf'  CK'I-VOM 

SliHfiOUTINfc  CNFV'<f'N(A»  FA*  NA.  X*  FX,  liX*  ':hT.  '■A<~F) 

C  I^TE^<■FOLATI()^'  OF  F(X)  FOR  AW  X  ARRAY  USINO  AITaFW"S  RFTh()D  FOR 
C  txTwAHOLATinw  OR  INTERPOLATION  NFAR  THE  ENDS  OK  THK  GIVEl.  ARRAY  AN 
C  NtVlLLF“'<  N'F'ThOD  for  ALL  OTHF  R  INTERPOLATION. 

C  INPUT 

L  A  =  THE  ARRAY  OF  OIVFN  AHSCISSA  IN  ASCFNOINO  ORUFP. 

C  K  A  =  THF  array  of  MVEN  VALLiFS  OF  F(A).  FA  IS  COPPLFX. 

c  NiA  =  THE  NUN-hER  of  elements  In  each  OE  ThE  APRAyS  a  and  F(A). 

C  X  =  THE  values  of  X  IN  ASCENDI'lO  ORDER  FOR  rhICh  OF'E  WISHES  TO 
C  determine  F(X). 

C  NX  =  THE  NUMPER  OF  EtFMFNTS  IN  EACH  OF  THE  APRAYS  X  AND  FIX). 

L  tiPT  =  THF  odd  number  of  points  USED  FOP  INTERPOLATION.  NPT  CAN 

C  NOT  PE  GPF.ATFR  than  NA. 

C  KASE  =  I,  THE  PROGRAM  PROCEEDS  NORMALLY 

C  KASE  =  ?,  FIX)  IS  calculated  ONLY  FOR  THOSE  X  VALUES  THAT  REC.UIR 
C  AINA)  FOP  INTERPOLATION 
C  OUTPUT 

L  EX  =  THE  ARRAY  OF  VALUES  APPROXIMATED  FOR  FIX).  Fx  IS  COMPLEX. 
COMPLEX  FA,  Fx.  FUNCT,  POLY 

DIMENSION  AINA),  FAINA),  X(NX),  FXINX),  FUNCTIIS),  AH5CI1S),  DIFIl 
IS),  POLVIlS) 

V600  format  IIhO/SOH  THF  X  VALUES  ARE  NOT  ARRANGED  IN  ASCENDING  ORDER./ 
I/5X.4HI  =  f20.R,SX.7hxiI)  =  E20.<J,SX,ahJ  =  E20 .9, Sx , 7hA  I J)  =  EP0.9 
2///1‘*X,HX.20X,4HEIX>  ) 

RbQl  FORMAT  (SX,2E?0,9) 

960?  FORMAT  (IHO*  THERE  ARE  NOT  ENOUGH  POINTS  IN  THE  GIVEN  ARRAY,*) 

LOOP  =  I 

IE  INPT  -  IS)  3.  3.  9 
9  NPT  =  IS 

3  IF  INPT  -  NA)  P,  H,  4 

4  NPT  =  f'RT  -  2 

IF  INPT  -  1)  5,  5.  6 

5  PPINT  9602 

call  Exit 

6  IF  ILOOP  .FO.  2)  GO  TO  3 

7  LOOP  =  2 
GO  TO  3 

B  NPT?  =  NPT/2 

IF  IKASF  .FO.  2)  GO  TO  1? 

11  NSTART  =  1 
MX  =  1 

GO  TO  16 

12  NSTART  =  MA  -  3 

TEST  =  AINSTAPT)  ♦  (A(r, START  ♦  1)  -  AINSTART))/? 

IF  INX  .LT.  1)  GO  TO  170 

DO  14  I  =  1,  NX 

IF  IX(I)  -  TEST)  14,  14,  13 

13  MX  =  I 
GO  TO  16 

14  CONTINUE 
170  CONTINUE 

16  NSTOR  =  NA  -  1 

IF  INX  .lT.  MX)  (,(,  TO  17S 
UO  12E  1  =  MX.  NX 
IF  IX  (  I  )  -AID)  135.  IS,  10 
10  IF  IX(I)  -  AINA))  as,  20,  130 
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15  FX(I)  =  FA<1) 

GO  TO  125 
20  FX(I)  =  FA(NA) 

Go  TO  125 

25  IF  (NSTOP  .LT.  NSTAPT)  GO  TO  180 
DO  85  J  =  MSTAPT,  NSTOP 
IF  (X(I)  -  A(J) )  32,  3S,  30 
30  IF  (X(I)  -  A(J  ♦  D)  AO,  H5 

32  11  =  I  -  1 

PRINT  9600,  I,  Xd).  J.  A(J) 

IF  (II)  34,  34,  33 

33  PRINT  9601,  (X(N),  FX(N),  N  =  1,  11) 

34  CALL  FXIT 

35  FX(I)  =  FA(J) 

NSTART  =  J 
Go  TO  125 

40  FX(I)  =  FA{J  ♦  1) 

NSTABT  a  J  ♦  1 
60  TO  125 
45  NSTART  =  J 

IF  (ABS(X(I)  -  A(J))  -  ABS(X(I)  -  A(J  ♦  1)))  50,  50,  55 
50  JJ  a  J 
GO  TO  60 
55  JJ  a  J  ♦  1 
GO  TO  60 
85  CONTINUE 
180  CONTINUE 

60  IF  <JJ  -  NPT2)  135,  135,  70 
70  IF  (JJ  ♦  NPT2  -  NA)  80,  80,  130 

80  KK  a  JJ  -  NPT2  -  1 

90  IF  (NPT  ,LT,  1)  GO  TO  185 
DO  95  K  a  i,  nPT 
KK  *  KK  ♦  1 
FuNCT(K)  a  FA(KK) 

AeSC(K)  a  A(KK) 

95  DIF(K)  =  ABSC(K)  -  X(I) 

185  CONTINUE 

NTOP  =  NPT  -  1 
LL  a  1 

100  IF  (NTOP  ,LT,  1)  GO  TO  190 
DO  105  L  a  1,  NTOP 
LLL  a  L  ♦  LL 

105  POLY(L)  a  (FUNCT(L)*OIF(LLL)  -  FUNCT  (L*  1 )  *DIF  (L )  )  /  ( AHSC  (LI.  L  ) 
1  -  ABSC(LI) 

190  CONTINUE 

IF  (NTOP  -  1)  l^0,  120,  110 
110  00  115  M  a  1,  NTOP 

115  FUNCT(N)  a  POLV(M) 

NTOP  a  NTOP  -  1 
LL  a  LL  ♦  1 
GO  TO  lOO 
130  INC  a  -  1 

KK  a  NA  ♦  1 
GO  TO  140 
135  INC  a  1 
KK  a  0 

140  IF  (NPT  .LT,  1)  GO  TO  215 


1^0  l4S  ^  -  li 
=  Kl<  ♦  INC 
FlINCT  (►  )  =  f  A  (H^  ) 

AFSC  (l<  )  =  A  ((<  K,  1 

lAF  UIMK  )  =  AHSC  (K  )  -  A  (  I  ) 

a  IS  contimif 

NTOP  =  ’''PT  -  1 
LL  =  1 

ISO  IF  (NTOP  .LT.  1)  or.  TO  ,-^0 
no  ISS  L  =  1.  NT(JP 
LLL  =  L  ♦  LL 

IsS  H()LV(L)  =  (Fii^rr  ( 1 )  «r)iF  (LLi. )  -  funct  (l*  l )  *f)lF  (ll)  )  /  (  ahsc  (lll  ) 
1  -  AbSC (LL> ) 

aao  CONTINUE 

IF  (NTOP  -  1)  lan,  l?0,  160 
160  00  16S  M  =  1,  NTOP 

165  FUNCT(P)  =  POLY(M) 

NTOP  =  NTOP  -  1 
LL  =  LL  ♦  1 
GO  TO  150 

120  FX(I)  =  POLY(l) 

125  CONTINUE 
175  CONTINUE 
WtTURN 
ENO 


*OfcCK  FLEAF 

COMPLEX  FUNCTION  FLEAF(taAvE«  H] ,  H2»  XDt  OFLTAP) 

C  FLAT  EARTH  Theory,  reference  3 

C'.MPLFX  TEMP,  0.  7.  72.  77.  HjJEhF.  wfRF7.  WERE.  7WE»F»  OELYAR 
HU  =  H?  -  Ml 

TEMP  =  (0.707106761?.  -  0 . 707 1 067b 1 2 ) »50HT ( . 5*WA vF ) 

XD2  =  SaHT(XOJ 

0  =  -tEMP*HO/XD2 

7  =  tFmp»0ELTAP*XD2  .  O 

Z/  =  -  7 

ZI  =  AIMAG(Z7) 

IE  (ZI.LT.O.  .0«.  (ARS(RtAL (771 ) .LT.6.  .ANO.  7I.LT.6.))  Gn  TO  10 
72  =  77**2 

HWERF  =  (72  -  2.)/(77*(72  -  3.5)) 

GO  TO  12 

C  WERE  =  COMPLEMFNTARY  EPpOP  FUNCTION 
10  WERE?  =  WERE (77) 

HwEPF  =  77  -  0.5*WFRF77  (77*WERF7  »  (O..  -  0  .Sf,4lAR58)  ) 

12  7WEWF  =  7  ♦  HWERE 

FlEAF  =  (C0*7WFPF  -  O.S)7(7*7v.EWF  -  0.5) 

RETURN 

END 
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COMPLFX  function  WFRF(777) 

C  COPPLFMENTAPT  frpop  function 

Complex  z*  zzz*  zvt  v*  Z2*  c*  w*  s 

UlMtNSION  C(l?>f  W(S,  4) 
tOUlVAL FNCt  (S.  C(l?) , 

DATA  (C(l)  =  (.0.  -  .5F)4l895ft3S)  ) 

UATA  W/(1..0.0)  . 

1  (3.b7H7944l  171‘*4?3E-01  .b.OTlbTTOSBAUPJTt-On  » 

2  (1  .831b63888P73418E;-02,3.40na62170660662F.-On  . 

3  ( 1 .23409fi040866786£-04,2.011S73170376004f-On  « 

4  (1.1?63'517471<>2646E-07.1.4S9S3S899001528E-01)  . 
b  ( 4. 27583576 1 bb8070E-0 1 .0.0000000 00 OOUO OOF *00)  . 

6  (3. 047.42n52569126E-01. 2. 0821893820283166-01) . 

7  (1.402395813662779E-01.2.222134401798991E-01)  . 

8  (6. 53177772a904697E-02. 1.7391831541634906-01)  ♦ 

9  (3.6281456489988646-02.1.3583895100065516-01 ) . 
A  ( 2. 553956763 105058E-0 1.0. OOOOOOOOOOOOOOOE.OO ) . 
b  (2. 184926152748907E-01. 9. 2997809392601866-02) « 
C  ( 1 .4795275951201586-01 .1 .31 1797170842178E-01 ) . 
0  (9. 27l07664264i.332E-02. 1.2831696222826156-01) » 
6  (5. 9686929610445906-02. 1.1321005612448826-01) » 
6  (1.7900 115 11 81 39306-0 1.0. 0000000000000006 *00)  . 
G  (1.64261 13639298616-01.5.0197135135249666-02) . 
H  ( 1 .3075746966985226-01 .8.1 112650477454726-02) . 
I  (9. 64C2505583044396-02. 9. 1236326004212586-02)  . 
J  ( 6. 9790 96 16496‘.750F -02. 8.9340000240 3646 16-02)/ 

XX  s  oEAl (7ZZ) 

YY  S  AIMAG(Z77) 

X  =  ABS(XX) 

V  =  A8<;(YY) 

2  -  CMPLXIX.  Y) 

L72  =  0 

16  (X  .G6.  4.5  .OP.  Y  .(>6.  3.5)  GO  TO  100 
1  =  X  .  .5 
J  =  Y  .  .5 

V  =  CMPLX (FLOAT (  I  )  ,  FL0AT{J)> 

ZV  =  7  -  V 

Cl2)  i  -(I  ♦  1.  J  ♦  1) 

Al  =  n. 

00  10  I  =  3.  12 
AI  =  AI  -  .5 

C  ( I )  =  (  Y»c ( I  -  1 )  ♦  C( I  -  2) )/AI 

10  CONT1n(jE 
J  =  12 

00  11  1  =  2.  11 

J  =  J  -  1 

11  S  =  S*7V  ♦  C(J) 

20  IF  (YY  .06.  0.)  GO  Tn  30 
IF  (L72  .60.  0)  72  =  7»7 
S  =  2.*CEXP(-72)  -  5 
IF  (XX  .OT.  0.)  S  =  CONJG(S) 

GO  TO  200 

30  IF  (XX  .LT.  0.)  5  =  CONJG(S) 

200  WtwF  =  S 
RETURN 
100  L72  =  1 


48 


7?  =  7*7 

S  a  7*(  (0.«0.461313S279)/(/2  -  0.  l9C)lti3bO'?2)  ♦  (O,, 
1/(72  -  1 . 7«4<t927<*8fij  *  (0..0.002«ft38938/<*8) /(72  -  ' 
00  TO  20 
END 


•deck  gang 

Function  cang(7) 

C  PUKPOSEi 

C  COMPUTES  THE  ARGUMENT  OF  A  COMPLEX  NUM^FP  7  SUCH  THAT 
C  -PI  .LT,  THETA  .LF.  PI. 

COMMON  /PIS/TWOPI  .PI.HAFPI.'lPTPI 

COMPLEX  2 

X«REAL(?) 

VsAIMAGlZ) 

IF  (X)  20*30*10 

10  CAN6sATAN2(r*X) 

RETURN 

20  IF  (Y.NE.O^)  00  TO  10 

CANGsPI 
RETURN 

30  IF  (Y.GT.O.)  GO  TO  40 
IF  (T.LT. 0.1  GO  TO  so 
CANG=0. 

RETURN 

40  CANGsHAFPI 
RETURN 

so  CANGs-HAFPI 
RETURN 
END 


0.09V9P21614H) 
-.52534374379) ) 


49 


•DECK  INDF 

SUBROUTINE  INDF  (HM,XX»E0RH,F ) 

C  THIS  SUBROUTINE  CALCULATES  THE  INDUCTION  FIELDS  FOR  E  SUB  R  AND 
C  H  SUB  PHI 

C  THESE  INDUCTION  FIELDS  APE  FOR  POSITIVE  TIME  FUNCTION 
C  EORH  =  1.  GIVES  INDUCTION  FIELD  FOR  E  SUB  R  (TOWEL  BAR) 

C  EORH  =  0.  GIVES  INDUCTION  FIELD  FOR  H  SUB  PHl  (LOOP) 

C  FZ  =  INDUCTION  FIELD  FOP  E  SUB  R 

C  FH  »  INDUCTION  FIFLD  FOR  H  SUB  PHI 

COMMON  /ZOTA/ARPAY ( 15) 

COMMON  /PIS/TWOPI *RI tHAFRI ,OPTRI 
DOUBLE  THETA0.SINTH,C0STH,R,C0NS 
COMPLEX  FZ.Fh.F 

A  =ARRAY(3)  5C  =  AR«AY(9)*l.E*9  SwAvE  =  ARRAY(IO) 

IFCXX.LE.O.  )3.4 

3  PRINT  9500 

9500  FORMAT  (//*  IN  INDF,  DISTANCE  IS  ZERO  OR  NEGATIVE,  XX  »  *,F20.10) 
CALL  EXIT 

4  THeTAO=l.O0*XX/A 
SINTH=OSIN(THfTAO) 

COSTH=OCOS(THFTAD) 

R=A*HH 

CONS=A*R*SINTH**2 

02=r*r*a*a-<;,*a*r*costm 
IF(02.GT.O,)  60  TO  30 
00=XX 
02sXX*XX 
GO  TO  40 
30  CONTINUE 
00*S<JRT<D2) 

40  continue 

03=00*n2 

04*02*02 

F2R*-2.*C0STH/D0*3.*C0NS/n3 

FZI= ( WAVE*C0NS*2.*C0STH/WAVE ) /02-3.*C0N5/ <04*WAVE ) 
F2=CMPLX(FZR,F7I) 

FZ«FZ/CMPLX (0,,WAVE) 

FHR*R*S INTH*00/03 

FHI*WAvE*R*SINTM*D0/D2 

FH*CMPLX(FHI,-FmR) 

Cla2.E-7*TwOPI*WAVE*C 

FhsFH/CI 

IFlEORH) 1 ,? 

1  F*FZ 
RETURN 

2  FsFH 
RETURN 
END 
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\  I 

^  MISSION  h 

\  of  V 

^  Rome  Air  Development  Center  C 

^  RAt?C  ptam  and  execota*  azAzoAch,  dzvzZopm&nt,  tzii  and  £ 

»  ^ztzcXzd.  aaqali,AM.on  pKogaam  -Ln  iuppoaX  0){  Command,  ConXAol  Q. 
g  CommyU.cationi>  and  JntzlUgzncz  (C^j)  acXlviUzA.  Tzchnical  ? 
i  and  ang-tneeAxng  iuppoAt  within  oavu  tzchnical  compztencz  C 
S  i^  pAovidzd  to  ESV  PAogaam  O^^iczi  (PO4}  and  otkzx  BSV  ^ 

^  edemen^.  Tkz  pai.nci.paZ  tzchnicaZ  mi&iion  oAeod  oat  \ 

1  cormanicatcom ,  eZectaomagnctic  guidance  and  contAoZ,  iua.-  A 

^  v&iZZancc  o^  gaound  and  aeAo^paez  obizeti,  intzltigzncz  data  ^ 
»  coZZzetion  and  handZing,  in^oamation  ^yitejn  tzeknoZogy,  v 

2  iono&phzAic  paopagation,  iotid  &taZz  ici-znezA,  micaowavz 
j;  phyiici  and  zZzetAonic  aztiabiJUty,  maintainabiZitu  and 
S  compatibiZity. 


Printtd  by 

Unitad  Statat  Air  Pore* 
Hantcom  AFB,  Mats.  01731 


